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ABSTRACT 

Most planetary systems are formed within stellar clusters, and these environ¬ 
ments can shape their properties. This paper considers scattering encounters be¬ 
tween solar systems and passing cluster members, and calculates the correspond¬ 
ing interaction cross sections. The target solar systems are generally assumed to 
have four giant planets, with a variety of starting states, including circular or¬ 
bits with the semimajor axes of our planets, a more compact configuration, an 
ultra-compact state with multiple mean motion resonances, and systems with 
massive planets. We then consider the effects of varying the cluster velocity dis¬ 
persion, the relative importance of binaries versus single stars, different stellar 
host masses, and finite starting eccentricities of the planetary orbits. For each 
state of the initial system, we perform an ensemble of numerical scattering ex¬ 
periments and determine the cross sections for eccentricity increase, inclination 
angle increase, planet ejection, and capture. This paper reports results from over 
2 million individual scattering simulations. Using supporting analytic considera¬ 
tions, and fitting functions to the numerical results, we find a universal formula 
that gives the cross sections as a function of stellar host mass, cluster velocity 
dispersion, starting planetary orbital radius, and final eccentricity. The resulting 
cross sections can be used in a wide variety of applications. As one example, we 
revisit constraints on the birth aggregate of our Solar System due to dynamical 
scattering and find fV< 10^ (consistent with previous estimates). 

Key words: planets and satellites: dynamical evolution and stability - plane¬ 
tary systems 
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A large frac tion of planetary systems form within stel¬ 
lar clusters (iLada fc Ladal I 2003 I : iPorras et al.l 120031 ') and 
these birth environments can in fluence their resulting 
prop erties (e.g., see the reviews of lAdani 3 l2010l:IPfaIzn^ 
120131 '). One potentially important process occurs when 
binary systems — and single stars — fly past solar sys¬ 
tems and disrupt the orbits of their constituent planets. 
This t ype of scattering interacti on has been studied in the 
held i Laughlin Ada msI 2000ll. and within young em¬ 


bedded clusters (e.g., AdamsetalJl200^ Malmbereet_^ 
200^j201lljBolevet_^ 201^ Duhes_&Kruruholz 20]A 


Chatteriee et al.r 20121 : Hao et ahT 2013l :~ Pacucci et al.l 

2013l l. where the latter results can be used to provide 
constraints on the possible birth environment of our own 
solar system ('e.g.. lAdams &: LaughlirJl200ll : [Hester et al.l 


derived for the birth aggegate of the solar system depend 
on many variables, including assumptions made about 
the cluster properties, any other constraints imposed on 
the problem, and the interaction cross sections. 

This present study focuses on the cross sections 
themselves, and expands previous work to include a much 
wider range of parameter space; the implications for the 
solar birth environment are then briefly considered at the 
end of the paper. For studies concerning our solar system, 
most previous work has calculated the cross sections for 
this mode of disruption by considering the initial orbits 
of the giant planets to have their present-day values of 
semimajor axis. However, some recent work suggests that 
our solar system may have begun in a more compact con- 
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figuration dGomes et al.ll2005l : iTsiganis et al.ll2005l ). and 
the planets may not have reached their present-day or¬ 
bits until the solar system reached an age of hundreds of 
millions of years. One motivation for this present study is 
thus to determine cross sections for solar system disrup¬ 
tion for more compact configurations. Note that the sign 
of the effect is not obvious a priori-. The geometrical cross 
section of the compact solar system is smaller, and hence 
implies a smaller interaction cross section. However, the 
decreased relative separations of the planets allow for in¬ 
creased planet-planet interactions, which could result in 
more disruption from the passing stars; in addition, the 
close spacing in compact solar systems allows for orbit 
crossing to occur for smaller values of eccentricity. 

In some compact configurations of the solar system, 
the giant planets can be at or near mean motion res¬ 
onance. This possibility leads to interesting dynamics: 
Mean motion resonances can protect planetary systems 
from disruption, and could thus lead to greater stabil¬ 
ity and smaller interaction cross sections. On the other 
hand, the mean motion resonances themselves are more 
easily compromised than planetary orbits — the poten¬ 
tial energy corresponding to the resonance angle being in 
a bound state is much less than the gravitational poten¬ 
tial energy of the planetary orbit. An important related 
question is thus to find the cross sections for passing stars 
(including binaries) to disrupt mean motion resonances. 
Planetary systems with disrupted resonances will usually 
retain their planets in the near term, although they could 
be subject to orbit instabilities over longer spans of time. 

In addition to compact solar system architectures, 
this paper considers a wider range of parameter space 
than previous studies. Part of this expanded scope is pos¬ 
sible due to increased computational capabilities. This 
present study includes results from more than 2 million 
individual numerical experiments that simulate a solar 
system interacting with a passing binary (or single star). 
For each choice of solar system architecture and each 
choice of the background parameters for the encounters, 
we run a large ensemble of Me simulations (where Me 
= 80,000 for most cases, but can be larger). The varia¬ 
tions that we consider for the target solar systems include 
compact configurations (described above), more massive 
planets, nonzero initial orbital eccentricities, and a range 
of masses for the central stars. Regarding variations in the 
background environment, this paper considers two main 
issues: We determine the effects of varying the velocity 
dispersion of the cluster stars, and we compare the rela¬ 
tive sizes of the scattering cross sections for single stars 
versus binaries as they interact with planetary systems. 

This paper is organized as follows. We formulate our 
approach to calculating the interaction cross sections in 
Section O The resulting cross sections are then given in 
Section [S] which provides (cr) for increases in eccentric¬ 
ity, increases in the spread of inclination angles, planet 
ejection, planet capture, and changes in semimajor axes. 
Results are also presented for increasing orbital eccen¬ 
tricities up to orbit-crossing configurations and compares 
the efficacy of passing single and binary stars. These re¬ 
sults are given as a function of solar system architecture, 
velocity dispersion of the cluster, and mass of the host 
star. Over much of the parameter space of interest, the 


cross sections display a nearly self-similar form. Section 
[4] presents a scaling analysis that shows how the results 
scale with velocity disperion, stellar mass, and starting 
semimajor axis. As an application. Section [5] revisits the 
possibile dynamical constraints on the birth cluster of 
the solar system. In order to assess the level of disrup¬ 
tion, one also needs the rate of close encounters in young 
stellar clusters. These rates have alre ady been calculated 
for a wide range of cluster properties llAdams et al.ll200^ : 
IProszkow fc AdamsIboOfll 'l and are used herein. The pa¬ 
per concludes, in Section[6l with a summary of our results 
and a discussion of their implications. 


2 FORMULATION OF THE PROBLEM 

One useful way to specify the effects that passing stars 
can have on planetary systems is to define cross sections 
of interaction. For example, the scattering interactions 
could eject a planet, increase the eccentricity, change the 
semimajor axis, and/or perturb the inclination angle of 
the orbit. For a given type of disruption, a solar sys¬ 
tem presents an effective target area for being disrupted 
by passing stars. With this definition, the effective inter¬ 
action rate T for disruption is then given by the usual 
formula 

r = n*(cr)(u} , (1) 

where n* is the mean density of stars in the environment, 
{v) is the mean relative velocity between systems, and 
(cr) is the cross section for the given mode of disruption. 
We note that the background environment determines 
the stellar density n* and the distribution of relative ve¬ 
locities. As outlined below, the relative velocities follow a 
Maxwellian distribution characterized by the expectation 
value {v). The interaction cross section depends on this 
velocity distribution, so that we actually calculate the 
quantity (cr)„ = {av)/{v), where the subscript denotes 
that the cross section depends on the velocity expecta¬ 
tion values (v). For ease of notation, however, we drop the 
subscript for the remainder of the paper. In young em¬ 
bedded clusters, we expect n* ~ 100 pc“® and (v) ~ 1 — 2 
km/s; in the field (in the solar neighborhood) these quan¬ 
tities have typical values n* ~ 1 pc“® and (v) ~ 30 — 40 
km/s. Because of the velocity dependence of the cross sec¬ 
tions, solar systems in the field (with high fly-by speeds) 
are, on average, less affected by passing stars. 

To calculate the cross sections for interactions, we 
adopt the following approach. First we must specify the 
configuration of the solar system that will be targetted 
for disruption (for example, we can use the current set 
of four giant planets in our solar system, with their cur¬ 
rent masses and semimajor axes, all in orbit about a solar 
mass star). Next we must specify the background envi¬ 
ronment, which determines the distribution of relative ve¬ 
locities. For most of this work, we focus on the case where 
the target solar system encounters binaries. We then per¬ 
form a large ensemble of numerical simulations, where 
the input parameters are specified according to a Monte 
Carlo scheme. The results are then used to calculate the 
probability of various outcomes and the corresponding 
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cross sections (for further deta i l, see Lauehlin fc Adatni 
l200fll : lAdams fc LaughlirJl200ll : lAdams et al.ll2nn6h . 

In principle, the Monte Carlo sampling scheme 
should sample all possible encounters between binaries 
and the target solar system, including those with large 
impact parameters. In practice, however, only sufficiently 
close encounters have a non-negligible chance of affecting 
the planetary orbits. In order to conserve computer time, 
we thus make the following limitation. We treat the semi¬ 
major axes a of the binaries on a different footing than 
the others: The values of a are sampled uniformly out to 
Umax = 1000 AU (more than 30 times the size of Nep¬ 
tune’s orbit in our solar system). For a given value of a, 
we then limit the possible range of impact parameters to 
fall within an area given by Aq = B-koI . With this sam¬ 
pling scheme, the cross section of interaction, for a given 
type of disruption event, is given by 


r^tmax 

(a) = / p{a)da foia) 
Jo 


( 2 ) 


where p(a) is the probability distribution for binaries hav¬ 
ing a semimajor axis a. The factor foia) represents the 
fraction of all encounters (within the pre-determined area 
Ao = BtvC) that results in the outcome of interest. Note 
that the maximum allowed value of the impact parameter 
varies with a and is given by Wmax = CBa. 

The formulation of equation m can be understood 
as follows: Consider a given outcome of interest, say, the 
ejection of Neptune. We only consider fly-bys that take 
place within the area Ao = Bnal, where a is sampled uni¬ 
formly. If every encounter within this area leads to the 
ejection of Neptune, and all encounters outside this area 
(which are not computed) have no effect, then fo = 1; 
the probability factor p{a) corrects for the actual distri¬ 
bution of binary semimajor axis, and one can see that 
equation provides the correct effective cross section. 
In practice, of course, only a small fraction of encoun¬ 
ters lead to the ejection of Neptune so that fo 1. 
As long as we choose the factor B large enough, we are 
ignoring only distant encounters that have little contri¬ 
bution to the cross section. Nonetheless, since B is finite, 
this procedure leads to a lower limit on the cross sec¬ 
tion. We have run convergence tests with ever-increasing 
values of B and find that B = 100 is large enough to in¬ 
clude essentially all relevant encounters. In most of this 
work we thus use B — 100, which provides a good com¬ 
promise between computationa l speed and accuracy. For 
comparison, our previous work (iLauehlin fc Adamj2000l : 
lAdams fc Laughlinll200lh used the smaller value B — 4, 
so that the reported cross sections (again presented as 
lower limits) were smaller than those obtained here by 
a factor of ~ 2. This present treatment thus provides a 
more complete accounting for wide binaries and results 
in a greater lower bound on the true cross sections. 

The distribution p{a) is determined by the observed 
binary period distribution, which is nearly uniform in the 
quantity log a, but has a broad peak centered at period 
P = 10^ days, which implies a ~ 42 AU for solar type 
stars touQuennov fc Mavorlll991^ . 

Within the scheme outlined above, encounters be¬ 
tween a given solar system and a passing binary are spec¬ 
ified by a large number of input parameters: We must 


specify the properties of the binary, including its semi¬ 
major axis a, orbital eccentricity es, the masses of the 
two stars Mi* and M2*, and finally the phase of the bi¬ 
nary orbit Oh at the start of the encounter. The orbital 
eleme nts (a, Cb) are sampled from their observed distribu¬ 
tions llDuauennov fc Mayor![l99 if) . Similarly, the stellar 
masses are sampled from a log-normal form of the stellar 
initial mass function (co nsis tent with that advocated by 
lAdams fc Fatuzz3ll996l and IChabrieillioO^ . Both mem¬ 
bers of the binary are sampled independently from the 
distribution and the stellar masses are limited to the 
range M* = 0.07 — IOMq. As a result, we exclude brown 
dwarfs and the very largest stars (which are both rare 
and tend to reside at cluster centers). The phase angle 
6h of the orbit is sampled uniformly over [0, 27r]. Next we 
must specify the incoming velocity Uoo of the solar system 
with respect to the binary center of mass; this speed is 
sampled from a Maxwellian distribution with a velocity 
dispersion Vb that characterizes the background environ¬ 
ment (e.g., a cluster). The remaining variables are the 
three angles {9, -i/i, (j>) necessary to specify the direction 
and orientation of the encounter, and finally the impact 
parameter zu. The impact parameter is chosen randomly 
within a circle of radius 10a centered on the binary center 
of mass (corresponding to the choice B — 100 in equation 

m)- 

Using a Monte Carlo scheme to select the input pa¬ 
rameters according to the distributions described above, 
we carry out a large ensemble of scattering simulations. 
For most cases we find that the number of simulations 
Me = 80, 000 is large enough to provide good statistics. 
The outcomes of these numerical experiments are then 
used to compute the fraction fo of disruptive encounters 
for a given type of outcome. The resulting errors due to 
incomplete sampling are typically 5 percent or less, but 
can be larger for rare events (e.g,. for planet ejections, 
the sampling errors are ~ 10 percent). 

Each simulation is thus an N-body problem. For 
most cases, N = 7, where the target system consists 
of four giant planets orbiting a host star and interacts 
with a binary. The equations of motion are int egrated 
using a Bulirsch-Stoer method (IPress et al.lll986[l . which 
allows for rapid integrations and high accuracy. Because 
we are interested in the planetary orbits, which only con¬ 
tain a small fraction of the total energy of the N-body 
system, the simulations must conserve the total energy 
to high accuracy in order to determine the final orbital 
elements. For example, the energy contained in the or¬ 
bit of Neptune, the least bound planet, is typically 10"^ 
times smaller than the binding energy of a binary, or the 
initial gravitational potential energy between the binary 
and the solar system. In practice, our individual simula¬ 
tions have an accumulated error of only one part in 10®, 
so that orbital changes are safely resolved. 


3 RESULTS FOR THE CROSS SECTIONS 

Using the formulation described in the previous section, 
we have performed several large ensembles of numeri¬ 
cal scattering simulations. Unless stated otherwise, we 
consider the solar systems to have four giant planets 
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and to interact with passing binary stars. We then con¬ 
sider a number of different solar system architectures 
for the starting states, as outlined below (see Table 1). 
To obtain reasonable statistics within the Monte Carlo 
scheme, the number of individual numerical experiments 
for each solar system architecture must typically be of 
order Me ~ 80,000. This choice produces relative er¬ 
rors (due to incomplete sampling) of order 5 percent or 
smaller. 

In the first set of simulations, we consider the tar¬ 
get system to be an analog of our present-day solar sys¬ 
tem. In this case, we place the four giant planets in orbit 
about a solar mass star and give the planets their cur¬ 
rent masses and semimajor axes. The eccentricities are all 
set to zero, however, so that we can measure the eccen¬ 
tricity increases produced by the scattering encounters. 
From the results of these experiments, we compute the 
cross sections for orbital disruption of each of the four 
planets (as outlined in the previous section). The results 
are shown as the solid blue curves in Figure [U which 
also presents the cross sections for a more compact start¬ 
ing configuration (described below). The error bars (not 
shown) due to incomplete Monte Carlo sampling corre¬ 
spond to relative errors with a root-mean-square (RMS) 
value of ~4.4%. 

In Figure [T] and throughout this paper, the cross 
sections for increasing the eccentricity to e = 1 incorpo¬ 
rate all of the ways that the planet can be removed from 
its solar system. These channels include [1] actually in¬ 
creasing the eccentricity to e ^ 1, which includes both 
hyperbolic orbits and planetary orbits that intersect the 
host star, [2] ejection from the solar system by increasing 
the kinetic energy so that the orbit is unbound, and [3] 
capture by one of the (two) passing stars. These chan¬ 
nels are not mutually exclusive, but the simulations are 
stopped after one of these events takes place. However, 
these channels only include ejection processes that hap¬ 
pen during or immediately after the encounter (we denote 
these processes as prompt ejections). In other cases, the 
planets are scattered into high eccentricity orbits, so that 
the orbits cross each other. With these configurations, in 
the absence of resonance, the planets will eventually ex¬ 
perience close encounters, which in turn lead to ejections 
or collisions (we denote this process as delayed ejection). 
The cross sections for delayed ejections will be considered 
later. 

For comparison, we also present the results from a 
series of numerical experiments using a more compact 
orbital architecture (shown as the red dashed curves in 
Figure [T]), which is m otivated by the Ni ce model of so¬ 
lar system formation l| Gomes et al .1120^ 1. Although the 
Nice model has a number of variations, one feature is 
that the giant planets could have formed with a more 
compact configuration than that of the present day. For 
this case, we fix the orbit of Jupiter at aj = 5.2 AU, and 
then let each successive planet have a semimajor axis 
that is larger than the previous one by a factor of 5/3. 
This evenly-spaced solar system thus extends out to only 
24 AU. The results, shown in Figure [T1 indicate that the 
cross sections for the compact configuration are some¬ 
what smaller than those obtained with the current semi¬ 



Figure 1. Cross sections for eccentricity increase for the cur¬ 
rent solar system architecture and for a more compact con¬ 
figuration movitated by the Nice model. For the current solar 
system (solid blue curves), the four giant planets are started 
with their current semimajor axes and zero eccentricity. For 
the compact configuration (dashed red curves), the planets 
are started with semimajor axes having a fixed ratio aj^ijaj 
= 5/3, where Jupiter {j = 1) is started at its present loca¬ 
tion aj = 5.2 AU. For both sets of cross sections, the curves, 
from top to bottom, correspond to Jupiter (bottom), Saturn, 
Uranus, and Neptune (top). Since the orbits start with zero 
eccentricity, the eccentricity increase Ae = e, where e is the 
post encounter eccentricity. 


major axes. For this compact solar system, the RMS er¬ 
rors (not shown) due to incomplete sampling are ~4.6%. 

Next we consider an even more compact orbital con¬ 
figuration, again motivated by the Nice model, where the 
four giant planets are in mutual mean motion resonance 
(MMR). In this case, we choose the starting semima¬ 
jor axes to have values of a = 5.88 AU (Jupiter), 7.89 
AU (Saturn), 10.38 AU (Uranus), and 12.01 AU (Nep¬ 
tune). With these semimajor axes, Jupiter and Saturn 
are in a 3:2 MMR, Saturn and Uranus are in a 3:2 
MMR, while Uranus and Neptune are in a 5:4 MMR 
(fo r further discussion of this initial state, and others . 


see Baty gm_fc Brownll2010l : iNesvornv fc Morbidelliir2012l : 


lUi fc Batygin 2014h . Note that the semimajor axis ra¬ 
tios do not imply period ratios with exact integer val¬ 
ues (although they are close). All of the orbital elements 
must be chosen properly to put the system in mutual 
MMR, and this requirement displaces the period ratios 
somewhat. Nonetheless, the resonance angles of the sys¬ 
tem (for all three planet pairs) are librating in the initial 
state, as required for MMR. With this initial state, the 
solar system is much more compact than at the present 
epoch, and the cross sections for interactions are smaller. 
This trend is illustrated in Figure [21 which compares the 
cross sections with those obtained for solar systems with 
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Solar System Architectures 


Configuration 

Jupiter 

Saturn 

Uranus 

Neptune 

Standard 

e = 0 

e = 0 

e = 0 

e = 0 

Compact 

a = 5.20 AU 

a = 8.67 AU 

a = 14.4 AU 

a = 24.1 AU 

Resonant 

a = 5.88 AU 

a = 7.89 AU 

a = 10.38 AU 

a = 12.01 AU 

Eccentric ^ 1 

e = 0.049 

e = 0.057 

e = 0.045 

e = 0.011 

Eccentric # 2 

e = 0.10 

e = 0.10 

e = 0.10 

e = 0.10 

Massive 

mp = Imj 

mp = Imj 

mp = Imj 

mp = Imj 


Table 1. Summary of solar system configurations. In the standard configuration (first line), the 
planets have the same masses and semimajor axes as those in our solar system but start with zero 
eccentricity. For the other configurations, the table entries list the initial values of the parameters 
that are different from those of the standard configuration (so that all of the unlisted parameters 
have their standard values). 



Ae 


Figure 2. Cross sections for eccentricity increase for the cur¬ 
rent solar system architecture and for a resonant configuration 
movitated by the Nice model. For the current solar system 
(solid blue curves), the four giant planets are started with 
their current semimajor axes and zero eccentricity. For the 
resonant configuration (dashed red curves), the planets are 
started with semimajor axes a = 5.88, 7.89, 10.38, and 12.01 
AU (for the analogs of Jupiter to Neptune). For both sets of 
cross sections, the curves, from top to bottom, correspond to 
Jupiter (bottom), Saturn, Uranus, and Neptune (top). Since 
the orbits start with zero eccentricity, the eccentricity increase 
Ae = e, where e is the post encounter eccentricity. 


the standard starting configuration. To leading order, the 
smaller cross sections obtained for the resonant architec¬ 
ture are a direct consequence of the smaller geometrical 
size. However, closer inspection of the results suggests 
that the cross sections are larger than the smaller size 
would imply (see the analysis of the following section). 
For example, the cross sections for changing the eccen¬ 
tricity of Uranus and Neptune are comparable. In this 
compact state, planet-planet interactions can be impor¬ 


tant and act to increase the the cross sections of Uranus 
(and Saturn) beyond the values obtained for more widely 
separated orbits. 

In addition to changes in the orbital elements of the 
individual planets, as shown in Figure [3 scattering in¬ 
teractions can remove solar systems from their resonant 
states. The energy required to remove a planetary system 
from resonance is much less than that required to eject 
a planet, or even to substantially change its orbital ele¬ 
ments. To address this issue, we have run an additional 
series of numerical simulations to determine the fraction 
of systems that are removed from their initial resonant 
state due to passing binaries. As before, the ensemble size 
A/fi « 80,000, although the simulations take longer be¬ 
cause the resonance angles must be monitored for several 
libration times after the encounters. The result of this set 
of experiments is the cross section for removing the solar 
system from its initial resonant state, namely 

(crjres « (2, 280, 000 ± 20, 800) AU^ . (3) 

This cross section is about 20 times larger than that re¬ 
quired to eject Neptune from the solar system in its nor¬ 
mal state, and nearly 40 times larger than the cross sec¬ 
tion to eject Neptune from the compact, multi-resonant 
state. If the removal of the system from resonance re¬ 
sults in orbital instability over longer time intervals, then 
the multi-resonant state could be more sensitive to dis¬ 
ruption from passing stars than the standard solar sys¬ 
tem architecture. We have carried out 70 longer-term in¬ 
tegrations for post-encounter systems and find that all 
but one are stable on time sca l es of ^ 1 Myr. Other 
authors (iBatygin fc Brownll2010l : iNesvornv fc Morbidellil 
|2012| ) also find that multi-resonant states can be unsta¬ 
ble due to perturbations (generally due to a planetesimal 
disk), and can eject planets, but more follow-up integra¬ 
tions are required to assess the probability of significant 
instability. 

The results reported thus far have all been calcu¬ 
lated for cases where the velocity dispersion Ub = 1 
km/s, a typical value for a n embedded cluste r envi¬ 
ronment llUada fc Ladal [20031 : iPorras et ahl l2003l ). Now 
we generalize the treatment by considering the depen¬ 
dence of the cross section on the velocity dispersion of 
the background environment. As is well known, interac¬ 
tion cross sections for high speed encounters, such as in 


© 2014 RAS, MNRAS OOO.IT1I231 


















(AU^) ^ a (AU^) 


6 Li & Adams 







Ae Ae Ae Ae 


Figure 3. Cross sections for eccentricity increase for the cur¬ 
rent solar system architecture over a wide range of velocity 
dispersions in the background cluster. The four giant planets 
of the solar system are started with their current semimajor 
axes and zero eccentricitiy. Each panel shows the cross sec¬ 
tions to increase orbital eccentricity for Jupiter (upper left), 
Saturn (upper right), Uranus (lower left), and Neptune (lower 
right). The velocity dispersions fall in the range from 1 km/s 
(uppermost curves in each panel) to 16 km/s (lower curves), 
and are equally spaced logarithmically (by factors of V2). 


Figure 4. Cross sections for a range of masses of the host star. 
Each case uses analogs of the four giant planets of our solar 
system, where the planets start with the current semimajor 
axes and zero eccentricitiy. Each panel shows the cross sections 
to increase orbital eccentricity for the analog Jupiter (upper 
left), Saturn (upper right), Uranus (lower left), and Neptune 
(lower right). The four curves in each panel correspond to four 
stellar masses, M* = 0.25, 0.5, 1.0, and 2.0 Mq, from top to 
bottom. 


the field (jLaughlin Adain^ 2QOOl) . are much lower than 
those in clusters (lAdams et alJ 200^), and the velocity 
dependence is relatively steep ( Adams fc Sperg^ I200!tI : 
iDukes fc Krumholj|2012 '). To study this dependence, we 
consider ensembles of numerical simulations with differ¬ 
ent values of velocity dispersion Vh- More specifically, we 
consider solar system starting with the current value of 
semimajor axes, and Ub in the range from 1 km/s to 32 
km/s, varied by factors of y/2 (so they are evenly spaced 
in a logarithmic sense). For the low end of this range of 
Ub, we can use the usual number Me = 80,000 of tri¬ 
als in the ensemble for each value of Ub- For the larger 
values of Ub, however, the cross sections are lower, and 
disruptive events are rare, so that we need larger values 
of Me to obtain good statistics (we find that the choice 
Me ~ 200, 000 is usually large enough). 

The interaction cross sections produced by this study 
are shown in Figure [31 where each panel corresponds to 
the results for one of the giant planets. The cross sections 
are plotted as a function of the post-encounter eccentric¬ 
ity e, for each choice of velocity dispersion vt- Figure [3] 
shows that the cross sections are almost evenly spaced 
in a logarithmic sense, with the lowest (highest) velocity 
dispersions producing the largest (smallest) largest cross 
sections. This finding suggests that the cross sections — 
to leading order — display a power-law dependence on 
the velocity dispersion. This claim is verified in the fol¬ 
lowing section. 


Next we consider the effect of changing the mass M* 
of the host star. Figure [4] shows the cross sections for 
systems with the current solar system architecture and 
varying stellar masses, from M« = 0.25 — 2.0 Mq. For 
these numerical experiments, the solar systems are all 
started with four planets that have the same masses and 
semimajor axes of the giant planets of our solar system. 
These analogs are labeled as ‘Jupiter’ through ‘Neptune’, 
although the host star can have a mass that differs from 
the Sun. As expected, the cross sections shown in Fig¬ 
ure [4] decrease as the stellar masses increases. Unlike the 
case of varying the velocity dispersion, however, the cross 
sections, considered as a function of eccentricity increase, 
do not display as much self-similarity: The cross sections 
decrease more steeply with eccentricity as the mass of 
the host star increases. Nonetheless, for a given value of 
eccentricity increase, cross sections for the four planets 
(with their four values of a) all show the nearly same 
(power-law) scaling with stellar mass. 

Notice that changing the stellar mass is (in one 
sense) akin to changing the planetary masses, because 
the mass ratios are the most important variables. How¬ 
ever, this association is not an equivalence: The masses 
of the passing binaries also enter into the problem, and 
their mass distribution is kept invariant. In addition, 
if the masses of the planets are increased to the point 
where the planet-planet interactions play a role, then self¬ 
excitation of eccentricity can produce larger cross sec¬ 
tions. This issue is addressed below where we consider 
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Figure 5. Cross sections for the solar system planets and 
varying initial eccentricities of the planetary orbits. For all 
cases, the four giant planets of our solar system are started 
with their current semimajor axes. The solid blue curves show 
the results for zero initial eccentricitiy; the dashed red curves 
show the results where the planets start with their current 
orbital eccentricities (e = 0.049, 0.057, 0.045, and 0.011); the 
black dotted curves show the results where the starting orbits 
all have e = 0.10. Cross sections are given for Jupiter (bottom 
curves), Saturn, Uranus, and Neptune (top curves), all given 
as a function of the post-encounter value e of the eccentricity. 


solar systems with larger planets. We expect interactions 
to be important in the regime where the angular mo¬ 
mentum exchange time scale between planets is com¬ 
parable to the encounter timescale. The exchange time 
scale can be determined, but the calculation is different 
for widely separated planets where the secular approx¬ 
imation is val id and for the resonant cas e (for further 
discussion, see lBatygin fc Morbidellill2013l) . 

For the starting configurations used thus far, the ini¬ 
tial orbital eccentricities of the planets have been taken 
to be zero. Given this choice, the resulting cross sections 
represent the cross sections for increasing eccentricity 
(which cannot decrease from its initial value). However, 
for the related problem of single stars interacting with 
binaries, an important difference arises between start¬ 
ing states where the binary has zero eccentricity and 
states where the binar y has small but finite eccentricity 
(iHeggie fc Rasiolll996l ). One might worry that the cross 
sections calculated herein could be affected by introduc¬ 
ing small starting eccentricities for the planetary orbits. 
We have explored this possibility by using two additional 
starting configurations for the solar system. In one case, 
the planetary orbits are started with their currently ob¬ 
served eccentricities, e = 0.049, 0.057, 0.045, and 0.011 
for Jupiter, Saturn, Uranus, and Neptune, respectively. 
In the second case, the planetary orbits are all started 
with a larger value of eccentricity e = 0.10. The result¬ 


ing cross sections are shown in Figure (5] along with our 
previous results with zero starting eccentricity. As shown 
in the Figure, all of the cross sections converge to the 
same values as long as the final eccentricity is moder¬ 
ately larger than the starting values. The difference be¬ 
tween results obtained starting with zero eccentricity and 
those where the orbits have their current eccentricity is 
modest. For the larger starting values e = 0.10, the cross 
sections for reaching e = 0.10 are enormous of course, 
much larger than the limits of the plot (and hence are 
not shown). Even for this starting state, however, the 
cross sections have almost converged to their “natural” 
values for e > 0.20, except for the case of Uranus; for this 
planet, the cross sections for eccentricity excitation only 
converge for e ^ 0.35. 

The results illustrated in Figure indicate that the 
problem of solar systems interacting with passing binaries 
is somewhat differ ent than that of singl e stars interact¬ 
ing with binaries ilHeggie fc Rasicj Il996l ). Starting with 
zero eccentricities has a larger effect in the binary-single¬ 
star setting. One difference between the two cases is that 
of symmetry: For a single star passing by a binary with 
zero eccentricity, the incoming trajectory is the same as 
the outgoing trajectory provided that the encounter is 
distant (so that the binary orbit can be considered as a 
ring of mass); this symmetry cancels some of the forc¬ 
ing. However, this symmetry is absent for solar system 
scattering, even when the planetary orbits are circular. 
The binaries that impinge upon the solar systems are 
themselves eccentric, where e is drawn from the observed 
binary eccentricity distribution (which favors high e). In 
addition, the solar systems have four planets, with dif¬ 
ferent orbital phases, and this property also breaks the 
symmetry (albeit to a lesser degree). Another difference 
between the two scattering problems is that the cross 
sections of this paper are averaged over an ensemble of 
different binary properties and differen t encounter pa¬ 
rame ters. The binary scattering results (iHeggie fc Rasiol 
1 19961 ) show that the the difference between finite eccen¬ 
tricity and circular orbits is largest for distant encounters, 
but the effect (the change in eccentricity) is largest for 
close encounters (see their Figure 2). The cross sections 
of this paper include both regimes, but the cross section 
is dominated by the close encounters where the results 
for e = 0 and e 7 ^ 0 are similar. As a result, starting 
the planetary orbits with non-zero eccentricity has only 
a modest effect on the cross sections considered in this 
paper (provided that one considers post-encounter eccen¬ 
tricities sufficiently larger than the starting values). 

Next we consider the effects of planetary mass on the 
scattering cross sections. The results are shown in Figure 
[^for the usual Solar System parameters and for an analog 
solar system where all of the giant planets have the mass 
of Jupiter (mp = Imj). Both classes of systems start 
with the same semimajor axes (the present-day values in 
our system) and zero eccentricity. The figure shows that 
the cross sections for increasing the eccentricities of Nep¬ 
tune and Uranus are largely unaffected by the increase in 
planetary mass, but the cross sections for Jupiter and 
Saturn are somewhat larger. Note that the cross sec¬ 
tions are plotted only for eccentricity values e ^ 0 . 20 . 
Within such a massive planetary system, small eccentric- 
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Figure 6. Cross sections for eccentricity increase in systems 
where the giant planets all have mass mp = Imj (dashed 
red curves). The cross sections for the current solar system 
architecture are shown for comparison (solid blue curves). In 
both cases, the planets are started with the current semimajor 
axes of the giant planets in our Solar System and with zero 
eccentricity. Cross sections are shown for analogs of Jupiter 
(bottom curves), Saturn, Uranus, and Neptune (top curves). 


Figure 7. Cross sections for increasing the spread of inclina¬ 
tion angles of the planetary orbits. All of the giant planets are 
started in the same plane; the quantity Ai is the total range of 
inclination angles of the four orbits after the encounters. Cross 
sections are shown for a variety of velocity dispersions, from 
Ub = 1 km/s (top curve) to Ub = 16 km/s (bottom curve), 
where the values are evenly spaced logarithmically (by factors 

of V2). 


ities (e ~ 0.10) are easily excited by planet-planet in¬ 
teractions; as a result, the cross sections for eccentricity 
increase — as determined through our numerical scheme 
— are extremely large and are not plotted in the figure. 

The numerical results for the cross sections can be 
understood as follows: To leading order, we often expect 
the planets to act as test particles, so that the cross sec¬ 
tions should not be sensitive to the planetary masses. For 
sufficiently massive planets, however, an increase in the 
eccentricity of one planet can lead to significant pertur¬ 
bations acting on the other planets, thereby leading to 
increased eccentricity excitation. By increasing the mass 
of all of the planets to that of Jupiter, the resulting solar 
systems are more excitable. The largest increase in the 
cross sections, which occurs for Jupiter and for low eccen- 
tricties, is only a factor of ~ 2; most cross sections expe¬ 
rience smaller changes. These results are generally con¬ 
sistent with the idea that our Solar System is “full”, i.e., 
no additional planets and little additional mass can be 
added to the extant planets without rendering the system 
unstable. In fact, even the current s olar system is unsta¬ 
ble on sufficiently long time sca les (iBatvein fc LaughlirJ 
I 2 OO 8 I : iLaskar fc Gastineaull2009| ') . 

Another way in which planetary orbits can be altered 
by scattering encounters is by changing their inclination 
angles. For all of the simulations, we start the four gi¬ 
ant planets in the same plane (so that ij = is = iu = 
*iv = 0). After the encounters, the inclination angles of 
the four planets are, in general, nonzero. We define the 


post-encounter spread Ai of the inclination angles ac¬ 
cording to the expression 


Ai = max 



' Jfc 1 

JjJk J J ’ 


(4) 


where the Ji are the angular momentum vectors of the 
planetary orbits and where the indices run through all 
four of the giant planets. The resulting cross sections for 
increasing the spread of inclination angles is shown in 
Figure [71 The Figure shows the cross sections for a range 
of velocity dispersions of the background cluster, from Ub 
= 1 km/s to Ub = 16 km/s, where the values are spaced 
by factors of \/2. The cross sections are almost evenly 
spaced in the semi-logarithmic plot and have nearly the 
same shape as a function of Ai. These properties indicate 
that the cross section has a power-law dependence on Wb 
(see Sections]). 

In general, increases in the inclination angles are pos¬ 
itively correlated with increases in eccentricity. This re¬ 
sult is not unexpected, as changes in both orbital ele¬ 
ments correspond to disruption of the initial states. To 
illustrate this trend, in Figure [8] we plot the increases in 
the spread of inclination angle Ai versus the change in 
eccentricity (equivalently, the post-encounter eccentricity 
since Ae = e). The two variables are in fact well corre¬ 
lated, but the range of possible Ai values for a given 
eccentricity e = Ae is large. As a result, in the figure we 
plot the mean values of Ai averaged over a bin in Ae 
with a width of 5 = 0.05. The data show a well-defined 
correlation; for this choice of binning, the spread in the 
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Figure 8. Correlation between eccentricity increases and in¬ 
creases in the spread of inclination angles of the planetary 
orbits. All of the giant planets are started in the same plane 
with circular orbits; the quantity Af is the total range of in¬ 
clination angles of the four orbits after the encounters. Corre¬ 
lations are shown for the orbital elements changes of Jupiter 
(heavy dashed red curve), Saturn (black solid curve), Uranus 
(black dotted curve), and Neptune (heavy blue solid curve). 
For each planet, the inclination angle increases are binned over 
a range in Ae of width <5 = 0.05. Although the correlation is 
well-defined, the range of Ai for a given value of Ae is rel¬ 
atively large. The error bars (shown for the Neptune curve 
only) depict the standard deviations. 


inclination angles grows to about 80° as the eccentricity 
grows to unity. The four curves shown in Figure [8] corre¬ 
spond to the four giant planets. Note that the orbits of 
all four planets show the same general trend. 

The cross sections discussed thus far correspond to 
the immediate, post-encounter properties of the solar sys¬ 
tems. In addition to immediate ejection, however, the so¬ 
lar systems can be rendered sufficiently unstable so that 
they eject planets long after the scattering encounters are 
over. These longer term ejection events can be divided 
into (at least) two types. In the first — and most unsta¬ 
ble — case, the scattering encounters leave the planetary 
orbits with high enough eccentricity so that adjacent or¬ 
bits cross. Most orbiting-crossing systems will eventually 
eject one of their planets, provided that the system is not 
in a mean motion resonance; furthermore, perturbations 
due to stellar encounters are unlikely to place a plane¬ 
tary system in resonance. We address the effects of this 
type of instability by finding the cross sections for pro¬ 
ducing orbit-crossing planetary systems (see below). In 
the second case, systems with more modest eccentrici¬ 
ties can be unstable over long spans of time. In order 
to assess the effects of this latter class of outcomes, the 
post-encounter systems must be integrated over typical 
stellar ages (billions of years). This task is beyond the 



Figure 9. Cross sections for the ejection of at least one planet 
as a function of velocity dispersion Vh in the cluster. The target 
systems have four giant planets with the masses and semima¬ 
jor axes of our solar system bodies. Cross sections are shown 
for three cases: increases in eccentricity large enough to pro¬ 
duce orbit crossing (solid dark curve marked by green error 
bars), direct ejection of a planet (dotted red curve), either 
channel of ejection (solid blue curve). 

scope of this present work, but provides an interesting 
problem for the future. 

Using the results of our numerical experiments, we 
can calculate the cross sections for the scattering inter¬ 
actions to leave any two orbits with high enough ec¬ 
centricities to cross. For the case of the analog solar 
system, where the four giant planets have their cur¬ 
rent masses and semimajor axes, the resulting cross sec¬ 
tions are shown in Figure |9l Three sets of cross sections 
are shown as a function of the velocity dispersion Vh of 
the background cluster. The cross sections for the post¬ 
encounter system to have an orbit-crossing configuration 
are shown as the lower, green solid curve in the figure. For 
the calculation of this cross section, only systems where 
all of the planets are retained by the host star are in¬ 
cluded. The cross sections for the system to eject any 
planet (including those planets captured by the pass¬ 
ing stars) are shown as the red dotted curve. Finally, 
the total cross sections for ejection, including both direct 
ejection of a planet and/or crossing orbits, are shown as 
the blue solid curve in the figure. The error bars depict 
the uncertainties in the cross sections due to incomplete 
Monte Carlo sampling. Note that the cross sections for 
orbit crossing and the cross sections for direct ejection 
are roughly comparable, with the latter slightly larger 
(except at high velocity dispersion, where they are the 
same within the sampling uncertainties). The total cross 
section for ejection is thus larger than that for direct ejec¬ 
tion by a factor of ~ 2. This statement holds only for the 
current solar system architecture, but remains valid over 
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Figure 10. Cross sections for changes in the semimajor axis of 
the planetary orbits due to scattering encounters. The target 
systems are analogs of our Solar System, with the four giant 
planets initially in circular orbits with the current values of 
their semimajor axes. The plots shows the cross sections for 
relative changes (Aa)/a in the semimajor axis for the orbits of 
Jupiter (lower red curve), Saturn (green curve), Uranus (cyan 
curve), and Neptune (upper blue curve). 


the range of velocity dispersion shown here (tib = 1—16 
km/s). 

We can now compare the results for the standard 
solar system architecture with that of the more com¬ 
pact configuration motivated by the Nice model. Here 
we consider only the most compact version where the 
planets are in multiple mean motion resonances (see Fig¬ 
ure ij. The compact configuration is expected to have 
lower cross sections for direct ejection. But the orbits are 
closer together, so that less eccentricity excitation is re¬ 
quired to produce crossing orbits. On the other hand, 
the semimajor axes are smaller, which lowers the cross 
sections for eccentricity increase. We find here that the 
cross sections for orbit crossing are comparable, (a) = 
96, 500 ± 3750 AU^ for the standard configuration ver¬ 
sus (cr) = 92, 200 ± 3710 AU^ for the compact multi¬ 
resonant case. However, the cross section for direct ejec¬ 
tion is larger for the standard solar system by a factor of 
1.5, so that the total ejection cross section remains larger 
by a factor of ~ 1.25. 

Although the semimajor axes of planetary orbits are 
altered less dramatically than the eccentricities and in¬ 
clination angles during scattering encounters, the values 
of a are nonetheless affected. The possible variations are 
quantified in Figure 1101 which shows the cross sections 
for producing relative changes {Aa)/a in the semimajor 
axes of the four giant planets. This ensemble of numerical 
simulations uses the standard solar system architecture 
as initial conditions, where the planets have their current 
masses and semimajor axes. The velocity dispersion of 


the background cluster is taken to be Ub = 1 km/s. As 
expected, the cross sections are largest for Neptune (top 
blue curve) and smallest for Jupiter (bottom red curve). 
As a crude approximation, the cross sections are propor¬ 
tional to the starting semimajor axes of the planets (al¬ 
though closer inspection shows the scaling is somewhat 
less steep than linear). 

Scattering encounters can cause the semimajor axes 
to become either smaller or larger, corresponding to the 
loss or gain of orbital energy. However, Figure fTOl shows 
that the process is highly asymmetric, where the orbits 
are much more likely to become larger (gain energy) than 
to move inward (lose energy). The scattering encounters 
rarely reduce the semimajor axes by more than a factor 
of two. Moreover, the magnitude of the cross sections are 
relatively small. More specifically, the cross sections for 
changing the initial semimajor axes by 10% are roughly 
comparable to — but somewhat smaller than — the cross 
sections for ejecting a planet (compare Figures [T] and [TOl. 
One might think that cross sections for moderate changes 
Aa would be larger than those for ejection. However, the 
cross sections for changes in semimajor axis do not in¬ 
clude the ejections themselves, i.e., they are the cross 
sections for changing the semimajor axis with the planet 
remaining bound to its host star. For large changes in a, 
there is not much parameter space where a is increased 
but the planet remains bound (thereby leading to the 
values shown in Figure HSl)- Notice also that the figure 
does not show cross sections for overly small values of 
(Aa)/a; the cross sections become singular in the limit 
(An)/o —>■ 0, as marked by the vertical dashed line. 

The cross sections considered thus far correspond to 
interactions between solar systems and passing binaries. 
On the other hand, roughly half of the stellar popula¬ 
tion consists of single stars, so that the corresponding 
cross sections for singles must also be determined. Since 
we are primarily interested in a comparison between the 
cross sections for single stars and binaries, it is crucial to 
use the same sampling for all of the parameters in the 
problem. Toward this end, we use exactly the same pro¬ 
cedure as before (outlined in Section [^, but let the mass 
of the second star go to zero. In this limit, the other, sin¬ 
gle star automatically resides at the center of mass of the 
system (and the value of the binary eccentricity becomes 
irrelevant). The resulting cross sections for single stars 
interacting with solar system analogs are shown in Fig¬ 
ure As before, the initial solar systems consist of four 
giant planets with the masses and semimajor axes of the 
present day Solar System (but with zero starting eccen¬ 
tricity). Each panel shows the interaction cross sections 
for eccentricity increases for a given planet (as labeled). 
Results are shown for four values of the velocity disper¬ 
sion of the background cluster, i.e., Vh = 1, 2, 4, and 8 
km/s (ordered from top to bottom in each panel). 

Next we make a rough comparison of the cross sec¬ 
tions for single star interactions (Figure 1111) with those 
obtained earlier for binaries (e.g.. Figure [3|. The single 
star cross sections are smaller by more than a factor of 
two. Note that the binary systems are, on average, some¬ 
what wider than the size of the solar systems. As a re¬ 
sult, as a pair of stars passes by a solar system, it consists 
mostly of empty space but still provides (roughly) twice 
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Figure 11. Cross sections for eccentricity increase due to 
encounters with passing single stars. The target systems are 
analogs of our Solar System, with the four giant planets in cir¬ 
cular orbits with their current values of semimajor axis. Each 
panel shows the cross sections for a given planet, as labeled, 
where the curves correspond to varying velocity dispersions of 
the background cluster: v^, = 1, 2, 4, and 8 km/s (from top to 
bottom). 


the opportunity for interaction as a single star. One thus 
expects at least a factor of two reduction in the cross sec¬ 
tions for passing singles. The fact that the reduction is 
larger than a factor of two is thus significant and indi¬ 
cates that the dynamics of the binaries themselves must 
contribute. Further, as discussed in the following section, 
the cross sections for single stars exhibit a different de¬ 
pendence on the background velocity dispersion and a 
slightly steeper dependence on post-encounter eccentric¬ 
ity. 

For convenience, Table 2 collects the cross sections 
for the ejection and capture of all four planets. For each 
solar system configuration considered in this paper, the 
table lists two sets of cross sections, where the first line 
corresponds to planetary ejection and the second line cor¬ 
responds to planetary capture. The Standard Model (the 
first configuration in the table) represents the case where 
the four giant planets have the masses and semimajor 
axes of our current Solar System, the host star has mass 
M* = I.OM 0 , the velocity dispersion of the cluster Ub = 1 
km/s, and the interacting stars are binary. The first col¬ 
umn in the table labels the solar system configuration by 
the variable that differs from its standard value. The er¬ 
ror bars in the table are those due to incomplete Monte 
Carlo sampling. One way to assess statistical significance 
is through the ratio of the cross section to its sampling 
error. For the ejection cross sections, the mean value (av¬ 
eraged over the entire table) of this signal to noise ra¬ 
tio is ~ 14, so that the ejection cross sections are well- 
determined. Capture events are much more rare. For the 


capture cross sections, the mean value of the signal to 
noise ratio is only ~ 4. For the rarest events, captures 
with high cluster velocity dispersion, the cross sections 
are only defined at the factor of two level. 


4 ANALYSIS AND SCALING LAWS 

The cross sections found in the previous section display 
relatively simple dependences on the underlying variables 
of the problem: For example, for each type of solar sys¬ 
tem, the cross sections, when considered as functions of 
the post-encounter planetary eccentricity e, all display 
the same general shape. As a result, the functions (cr)(e) 
can (almost) be rescaled to find a universal functional 
form, where scaling factors take into account the initial 
semimajor axis a of the planet, the velocity dispersion 
Vb of the background environment, the stellar mass M*, 
and so on. The goal of this section is to understand the 
general scaling properties of the cross sections and to 
determine the extent to which they are self-similar. In 
general, self-similarity arises when physical scales are ei¬ 
ther mi ssing from a prob lem or do not contribute to the 
results llBarenblattll200ijj : we return to this issue at the 
end of the section. 

Even in the reduced case where we consider one 
planet at at time, the interactions considered in this pa¬ 
per involve four bodies (the host star, the planet, and 
two binary members). Unfortunately, four-body interac¬ 
tions are rather difficult to describe analytically to any 
reasonable degree of approximation. As a result, the goal 
of this section is relatively modest: Instead of building 
complicated analytical models for 4-body (and higher N- 
body) dynamics, we consider here basic physical princi¬ 
ples that can be used as motivation for scaling laws. We 
then combine these heuristic results with our detailed nu¬ 
merical determinations of the cross sections. The result 
is physically motivated htting formula that characterize 
the cross sections over the parameter space of interest 
{a,Vb,M,,e). 

To start the discussion, consider the simplest case 
where the the cross section for interactions is the geo¬ 
metrical cross section rra^ provided by a planet in its 
initial orbit. Further, we consider the planets to be inde¬ 
pendent of each other during the encounters. This cross 
section will be enhanced by gravitational focusing, so we 
can write down an heuristic expression for the cross sec¬ 
tion in the form 

(ct)0 OTTU^ (^1-I-, (5) 

where Vesc is the escape speed from the target system 
(at the location of the planet), and v^o is the asymptotic 
relative speed between the two systems. In order to pass 
within this cross sectional area, the interacting star (bi¬ 
nary) must be about the same distance from the planet 
as its host star, so that the planet has a chance of being 
ejected from its original solar system. This expression 
thus represents the escape cross section. The parameter 
a is a dimensionless constant of order unity and is in¬ 
cluded to encapsulate the uncertainties inherent in this 
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Cross Sections for Ejection and Captnre 


Configuration 

Jupiter 

Saturn 

Uranus 

Neptune 

Standard model 

15500 ± 1360 

34000 ± 2050 

72300 ± 3100 

113000 ± 3860 


812 ± 306 

2140 ± 531 

6040 ± 994 

11400 ± 1280 

Compact model 

18100 ± 1510 

32700 ± 2130 

57500 ± 2790 

93900 ± 3570 


915 ± 379 

2280 ± 607 

4380 ± 817 

11500 ± 1320 

Resonant model 

23900 ± 1810 

40200 ± 2440 

61100 ± 2990 

60100 ± 2900 


1240 ± 467 

2150 ± 569 

3430 ± 701 

3620 ± 738 

Massive planets 

24100 ± 1890 

38300 ± 2390 

77700 ± 3360 

105000 ± 3880 


1530 ± 579 

2170 ± 637 

4810 ± 932 

8480 ± 1190 

Ub = 2 km/s 

9170 ± 947 

14800 ± 1250 

29800 ± 1770 

45200 ± 2240 


391 ± 136 

635 ±173 

2600 ± 487 

6370 ± 903 

Ub = 4 km/s 

2980 ± 454 

6090 ± 776 

10600 ± 918 

15700 ± 1140 


270 ± 173 

607 ± 230 

1580 ± 486 

2430 ± 569 

i;b = 8 km/s 

1220 ± 258 

2580 ± 403 

4060 ± 506 

5830 ± 698 


130 ± 85 

134 ± 75 

239 ± 88 

624 ± 228 

i^b = 16 km/s 

181 ± 39 

607 ±113 

1220 ± 182 

2140 ± 252 


82 ± 52 

53 ± 39 

214 ± 83 

169 ± 69 

M. = 0.25Mq 

37400 ± 2195 

74800 ±3170 

138000 ± 4350 

196000 ± 5130 


766 ± 330 

3510 ± 742 

10300 ± 1240 

19500 ± 1720 

M. = O.SMq 

27600 ± 1830 

54000 ± 2630 

107000 ± 3740 

152000 ± 4460 


1730 ± 560 

3310 ± 755 

7470 ± 1060 

15100 ± 1490 

M« = 2.0Mq 

11000 ± 1170 

21700 ± 1720 

45700 ± 2420 

69300 ± 2980 


458 ± 234 

1590 ± 467 

4420 ± 844 

7410 ± 1110 

Single, Ub = 1 km/s 

3840 ± 651 

7100 ± 856 

17300 ± 1430 

30600 ± 1980 


135 ± 80 

587 ± 210 

2080 ± 480 

5090 ± 871 

Single, Ub = 2 km/s 

1620 ± 324 

3110 ± 429 

9030 ± 926 

13200 ± 1090 


168 ± 86 

236 ± 106 

1370 ± 423 

2810 ± 641 

Single, Ub = 4 km/s 

685 ± 177 

1300 ± 244 

3740 ± 531 

6790 ± 793 


116 ± 94 

117 ± 51 

360 ± 117 

1480 ±411 

Single, Ub = 8 km/s 

269 ± 103 

1090 ± 322 

1440 ± 286 

1880 ± 266 


21 ± 14 

23 ± 14 

157 ± 53 

374 ± 202 


Table 2. For each solar system configuration, as labeled in the left column, the top line lists the 
ejection cross sections and the second line lists the capture cross sections for each of the planets. 
The error bars due to incomplete Monte Carlo sampling are included. All cross sections are given 
in units of AU^. 


approximation. After inserting the expression for the es¬ 
cape speed, we obtain 

(cr)o = arra^ aula , (6) 

where we have replaced the asymptotic speed Voo with 
the velocity dispersion Vh of the clnster (or other back- 
gronnd stellar system) and we have defined the corre¬ 
sponding length scale I = GMt/v'^ (where I ~ 890 AU 
for Vb = 1 km/s). The final expression represents the 
limiting form, which is applicable when gravitational fo¬ 
cusing dominates, and implies a linear dependence of the 
cross section on a. Given this form, the cross section re¬ 
quires another length scale. In this problem, the orbit 
speed of the binary, the asymptotic speed Vaa of the en¬ 
counter, the orbit speed of the planet, and the velocity 
dispersion Ub are all roughly comparable (1-10 km/s). 
For example, the orbit of Neptune in our solar system has 
an escape speed of ~ 5.5 km/s, whereas the orbit speed 
of a binary at the peak of the period distribution is also 
~ 5 km/s. If we only have a single velocity V, then di¬ 
mensional analysis implies that the relevant length scale 


must hel = GM^/V'^, as given in equation ((6|; additional 
uncertainties can be absorbed into the dimensionless pa¬ 
rameter a. Finally we note that for a velocity dispersion 
Ub = 1 km/s, the gravitational focusing term dominates 
by a factor of 30. 

The limiting form of equation is linear in the 
starting semimajor axis a of the planet. To see how well 
this expression works, we plot the ejection cross sections 
of the planets versus semimajor axis in Figure fT^ As ex¬ 
pected, the ejection cross section is a nearly linear func¬ 
tion of the semimajor axis. This trend holds for solar sys¬ 
tems starting with the present-day semimajor axes (star 
symbols) and the more compact configuration where the 
semimajor axes are spaced by factors of 5/3 (open tri¬ 
angles). We also plot results for ultra-compact solar sys¬ 
tems in multiple mean motion resonance (open squares). 
In order to isolate the dependence of the cross sections 
on initial semimajor axis from planet-planet scattering 
effects, this latter case uses smaller planet masses (by a 
factor of 10), so they act more like test masses; to explore 
a wider range in a, we also take this compact system to be 
smaller by a factor of 1.35 compared to that considered 
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Figure 12. Cross sections for planetary escape versus the 
starting semimajor axis. The 12 points on the plot correspond 
to the four giant planets in each of three versions of the ini¬ 
tial solar system architecture. The symbols represent different 
starting states, including the semimajor axes of the present- 
day solar system (stars), a compact configuration with 5/3 
semimajor axis ratios (open triangles), and an ultra-compact 
solar system starting in multiple mean motion resonances 
(open squares). The red solid line shows the cross section in¬ 
dicated by the limiting form of equation the blue dashed 
curve shows the full form. 

in the previous section. The error bars delineate the un¬ 
certainty due to incomplete Monte Carlo sampling. Not 
only do the cross sections show nearly linear dependence 
on a, but the slope of the curve is predicted by the above 
analysis. The red solid (blue dashed) curve in Figure [12] 
shows the cross section predicted by equation for the 
limiting case (full form); for both cases, the characteristic 
length scale £ = 890 AU and the dimensionless parameter 
a = 7/50 

Next we consider the dependence of the cross sec¬ 
tions on the post-encounter eccentricity e (which is equiv¬ 
alent to Ae because the orbits start with zero eccentric¬ 
ity). For all four planets in all three types of solar system, 
the e-dependence is similar. Since the ejection cross sec¬ 
tions scale linearly with semimajor axis a (see Figure [T2ll . 
we scale the cross sections by dividing out one power of 
a. The resulting scaled cross sections are shown in Fig¬ 
ure m as a function of eccentricity e. In addition to the 
individual cases (shown as the light dotted curves), the 
average is shown as the heavy blue curve, where the error 
bars depict the standard deviation. This latter quantity 

^ In order to set the value for this dimensionless parameter, 
and others specified in this section, we generally search in in¬ 
crements of 10“^, find the value that gives the minimum RMS 
error, and then choose the nearest round number (ratio of rel¬ 
atively small integers). 



Ae 


Figure 13. Scaled cross sections versus eccentricity increase 
Ae (equivalently, the post-encounter eccentricity e) for the 
four giant planets in each of the starting architectures for the 
solar system. The individual cases are shown as light dotted 
curves. The heavy solid blue curve depicts the average, where 
the error bars depict the standard deviation. The straight red 
line shows the result for cross sections with a purely exponen¬ 
tial dependence on eccentricity. 

provides a measure of the spread in the values of the 
cross section over the various cases. The standard devi¬ 
ation varies from about 17% of the cross section at low 
eccentricity e = 0.10 to only about 9% at e = 1.0. 

The curves in Figure fT51 are nearly straight lines on 
the semi-logarithmic plot, so that the dependence of the 
cross sections on eccentricity is nearly exponential. For 
purposes of illustration, we use an exponential fitting 
function of the form 

= a7r.£exp[6(l — e)] , (7) 

where the first factor enforces consistency with the ejec¬ 
tion cross sections considered above. For the value b = 
4/3, we obtain a good fit to the calculated, scaled cross 
sections, as shown by the heavy red line in Figure fT^ Ex¬ 
cept for first point (e = 0.1), the exponential fit (straight 
red line) agrees with the average values (solid blue curve) 
to within about 3%, i.e., the difference is much less than 
the width of the distributions as measured by the stan¬ 
dard deviations. Another measure of the quality of the 
fit is provided the relative differences between the nu¬ 
merically determined cross sections used in constructing 
Figure [T3] and the exponential form given by equation 
0 ; the RMS of these relative errors is ~ 12%. 

Next we consider the effects of the velocity dispersion 
of the background cluster environment. As shown in Fig¬ 
ure [3] in the previous section, the cross sections vary with 
the post-encounter eccentricity with approximately the 
same functional form over a wide range of Ub- Only the 
leading cofhcient changes. Moreover, the uniform spacing 
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Figure 14. Scaled cross sections versus eccentricity increase 
Ae (equivalently, post-encounter eccentricity e) for a collec¬ 
tion of different velocity dispersions for the background. The 
starting state is taken to have four giant planets with the cur¬ 
rent semimajor axes. Cross sections are scaled by ja (see 
text). The individual cases are shown as light dotted curves, 
with include curves for each of the planets for un = 1 ^ 16 
km/s, equally spaced logarithmically (by factors of v^)- The 
heavy solid blue curve depicts the average, where the error bars 
depict the standard deviation. The red striaght line shows the 
result for cross sections with a purely exponential dependence 
on eccentricity. 

of the curves in Figure |3] indicates that the cross sections 
must have a power-law dependence on the velocity dis¬ 
persion Ub (to leading order). We have explored scalings 
with velocity dependence of the form (a) oc and find 
that the best fit occurs for 7 « 7/5. Using this choice 
of power-law index, we plot the scaled cross sections ver¬ 
sus post-encounter eccentricity in Figure 1141 where we 
include the linear a-dependence found previously (i.e.. 
Each light dotted curve in the figure shows 
the result for one planet and one choice of velocity dis¬ 
persion. The heavy blue curve shows the average over 
all of the curves, where the error bars depict one stan¬ 
dard deviation. The heavy straight red line represents the 
same exponential dependence given in equation 0 and 
used in Figure 1131 The RMS of the relative differences 
between the numerically determined cross sections and 
the curve given by equation 0 is ~ 13%. The cross sec¬ 
tion curves are thus self-similar to this level of accuracy. 
Furthermore, the dependence of the cross sections on ve¬ 
locity dispersion is nearly independent of the dependence 
on starting semimajor axis a of the planet. 

The dependence of the interaction cross sections on 
the mass of the host star is somewhat more complicated 
than for the other variables, as illustrated in Figure |4l 
As the mass M* of the star increases, the cross sections, 
considered as functions of eccentricity, become steeper. 



Figure 15. Scaled cross sections versus eccentricity increase 
Ae (equivalently, post-encounter eccentricity e) for solar sys¬ 
tems with different stellar masses. The starting state is taken 
to have four giant planets with the current masses and semi¬ 
major axes. Cross sections are scaled by ja (see text). 

The individual cases are shown as light dotted curves, with 
include curves for each of the four planets for four choices of 
stellar mass M* = 0.25 — 2.0 Mq (spaced by factors of 2). 
The heavy solid blue curve depicts the average, where the er¬ 
ror bars depict the standard deviation. The red striaght line 
shows the result for cross sections with a purely exponential 
dependence on eccentricity. 


The spacing of the curves in Figure |4] (for different stel¬ 
lar masses) grows with e, so that the curves are not self¬ 
similar. In spite of this complication, we can still ht the 
cross sections with a power-law function of stellar mass, 
although the accuracy of the approximation is not ex¬ 
pected to be as high as in the previous cases. We thus 
consider a scaling of the form (a) oc and vary the 

index /r to hnd the best fit. The choice p = 1/3 provides 
the lowest RMS of the relative error. Figure [15] shows 
the result by plotting the scaled cross sections ja 

(again including the linear dependence on semimajor axis 
a) as a function of post-encounter eccentricity. The light 
dotted lines show the individual (scaled) cross sections 
and the heavy blue curve shows the average. The error 
bars depict the corresponding standard deviation, which 
is larger than for the cases considered previously (com¬ 
pare Figure [ 15 ] with Figures [T3] and [T4ll . The heavy red 
staight line shows the same result as before (from equa¬ 
tion [Zl). The RMS error between the exponential line 
and the numerically determined cross sections is about 
20%. This larger error measure results from htting the 
cross sections with a power-law form, even though the 
results depart somewhat more from self-similarity. 

The cross sections for increasing the spread of incli¬ 
nation angles, considered over a range of velocity disper¬ 
sions, also show a nearly self-similar form (see Figure [7||. 
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Figure 16. Scaled cross sections for increasing the post¬ 
encounter spread Ai of the inclination angles of the planetary 
orbits. The starting states have the four giant planets orbiting 
in the same plane (Ai = 0). The cross sections are scaled by 
the velocity dispersion of the cluster with the relation . 

The individual cases are shown as light dotted curves. The 
heavy solid blue curve depicts the average, whereas the error 
bars depict the standard deviation. The heavy red curve shows 
the fitting function described in the text. 

This finding indicates that the cross section should scale 
with a nearly power-law dependence so that (cr) oc . 
Over the range Ub = 1 — 16 km/s, we find that the best 
fit occurs for 77 ~ 7/5. To illustrate how well this scaling 
law works, we plot the scaled cross sections as a 

function of sin(Ai) in Figure fTBl Each light dotted curve 
in the figure corresponds to the result of one choice of ve¬ 
locity dispersion. The heavy blue curve shows the average 
of the scaled cross sections, where the error bars depict 
the standard deviations. The mean size of the error bars 
corresponds to relative differences of ~ 6 %, so that the 
curves are self-similar to this degree of accuracy. Notice 
that the scaling exponent 77 « 7/5 for inclination angle in¬ 
creases as a function of velocity dispersion Ub is the same 
as the corresponding index for eccentricity increases. 

After the velocity dependence has been scaled out, 
the cross section for increasing the spread of inclination 
angles is a slowly varying monotonic function of Ai (see 
Figure [T6ll . If we consider x = sin Ai as the independent 
variable (instead of Ai itself), the cross section can be 
fit with an exponential function which is analogous to 
that used to describe the eccentricity dependence. More 
specifically, if we use the functional form 

{a)i = {a}o exp [bo (1 — sin Ai)] , (8) 

then the cross section for increasing Ai can be fit using 
the parameters b ~ 3/4 and ao ~ 166, 000 AU^. Note that 
the value of the index b used here somewhat smaller than 
that needed to fit the dependence of the cross sections 


on (post-encounter) eccentricity (compare with equation 
[ 7 ]). The fitting function from equation ([8]) is shown in 
Figure [16] as the solid red curve. The quality of the fit 
is reasonably good: The fitting curve falls within one 
standard deviation (marked by errorbars in the figure) 
of the mean for all of the range except the first point 
(x = sinAi = 0.1); alternately, the RMS of the relative 
error between the two curves is ~ 8%. However, the mean 
of the numerical results (blue curve) shows more curva¬ 
ture than the exponential fit (red curve), especially at 
small values of x. 

Although we could find a more complicated fitting 
function that has smaller RMS relative error, we use 
equation (|8| in order to compare changes in the spread of 
inclination angle with changes orbital eccentricity. If we 
equate the variable x = sin Ai with e, then equations © 
and (|8|) have the same general form. We can then compare 
the leading coefficients, which have values (cr)o ~ 166,000 
AU^ for Ai-dependence and aTvia ~ 120,000 AU^ for e- 
dependence, where we have used a = 30 AU to evaluate 
the latter expression. The cross sections for eccentricity 
increase and spread of the inclination angles thus display 
similar behavior. The leading coefficients agree to within 
~ 28% and we can make the following inexact analogy: 
An increase in Neptune’s eccentricity of Ae = 0.1 corre¬ 
sponds to changing the spread of the inclination angles 
(of all four planets) so that sinAi increases by 0.1. We 
can also make a rough association between increasing the 
spread of inclination angles to Ai ^ 90 degrees and the 
ejection of a planet (e > 1). Both of these events have 
(approximately) the same cross section and both involve 
order unity changes to the angular momenta of the plan¬ 
etary orbits. In addition, the cross sections for inclination 
angle increase and eccentricity increase scale with the ve¬ 
locity dispersion in the same manner (oc u^*). 

The association between changes in the variables 
sinAi and e provides an intriguing topic for additional 
work. To leading order, the canonical actions written in 
terms of the orbital elements have the forms 

F oc and Z oc sin^(i/2) . (9) 

The apparent relation between the two variables (as 
observed in the simulation results) could thus be evi- 
dence of an equipartition-like m ixing of the actions (see 
iLichtenberg fc Lieberm^ Il992f l. Although beyond the 
scope of the present paper, this issue should be explored 
further. 

We can extract a potentially important cross section 
from these results. The scattering interactions considered 
here can readily increase the spread of inclination angles 
of outer bodies in a solar system. On the other hand, the 
scattering events themselves have little effect on planets 
in tight orbits, such as th e multi-planet system s observed 
by the Kepler mission llBatalha et al.l [2013ll . However, 
the bodies in the outer solar system can have important 
long-term effects on the inner bodies provided that they 
are scattered into orbits with sufficiently high inclination 
angles. More specifically, if the inclination angles of the 
outer orbits are larger than 39.2°, th en the Kozai effect 
can operate llKozaill 19621 : lLidovlll962ll . and the inner por¬ 
tions of the solar system can be excited over the age of 
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Figure 17. Scaled cross sections versus eccentricity increase 
Ae (equivalently, post-encounter eccentricity e) for solar sys¬ 
tems interacting with single stars. The starting states have 
four giant planets with the current masses and semimajor axes 
of our Solar System. Cross sections are scaled by the factor 
/a. The individual cases are shown as light dotted curves, 
which include the four giant planets and four values of veloc¬ 
ity dispersion of the background cluster: v\^ = 1, 2, 4, and 8 
km/s. The heavy solid blue curve depicts the average, where 
the error bars depict the standard deviation. The red striaght 
line shows the result for cross sections with an exponential 
dependence on eccentricity. 


the systems. Combining this requirement with the results 
of our numerical simulations, we find that the cross sec¬ 
tion for scattering a solar system into a state where the 
Kozai effect can operate is given by 

™ (srai) (uSt:) 

where flout is the semimajor axis of the outermost planet 
of the system. Note that the requirement of large mutual 
inclination is necessary but not sufficient for the Kozai 
effect to play a role. The Kozai effect is a highly fragile 
type of interaction because it involves libration of the ar¬ 
gument of periastron, and this quantity can be subject to 
ma ny other sources of precession (for further discussion, 
see lBatygin et aDl201ll b We also note that this form for 
the cross section (equation [TO]) involves some extrapo¬ 
lation: The numerical simulations were carried out pri¬ 
marily for the architecture of the current solar system. 
Nonetheless, the outermost planet is always the most af¬ 
fected by fly-by interactions, and the cross sections scale 
linearly with semimajor axis to a good approximation. 

Next we consider the scaling behavior of the cross 
sections for interactions with passing single stars. As for 
the case of binary systems, we expect the cross sections to 
scale nearly linearly with the semimajor axis a of a given 
planet. In addition, the nearly equal spacing on the loga¬ 


rithmic plot of Figure [Tl] indicates that the cross sections 
should display power-law dependence on the velocity dis¬ 
persion, such that (cr) oc . The velocity dependence 
for these single star cross sections is moderately less steep 
than those found earlier for binaries; the optimal value of 
the index 7 s ~ 6/5, which is somewhat smaller than the 
value for binary cross sections 7 « 7/5. After scaling out 
the semimajor axis and velocity dispersion, the reduced 
cross sections are shown in Figure 1171 The light dotted 
curves show the scaled values for given planets and values 
of Vh (which lie in the range 1-8 km/s). The heavy blue 
curve shows the mean over the entire collection and the 
error bars denote the standard deviations. These error 
bars correspond to an average relative error of ~ 15%, 
which is comparable to, but somewhat larger than that 
found for the binary cross sections. 

The scaled cross sections shown in Figure [T7] for 
single star interactions show a nearly exponential de¬ 
pendence on the post-encounter eccentricity. Although 
this behavior is analogous to that found for the binary 
cross sections, the slope of the exponential is somewhat 
steeper. Here we consider a fitting function of the form 

(cr)si„gle = (cr)o j [^=(1 - e)] ’ (11) 

where we obtain a good fit for (< 7)0 = 1000 AU^ and 
ba = 8/5. The resulting fit is shown as the red straight line 
in Figure [TtI The RMS difference between the expression 
of equation 113 and the numerically determined, scaled 
cross sections for single stars is only ~ 8%. 

Now we can compare the cross sections for single 
stars with those for binaries. The comparison is compli¬ 
cated by the different scalings of the two cases with veloc¬ 
ity dispersion and the different exponential laws for the 
eccentricity dependence. To fix ideas, consider the bench¬ 
mark case where the velocity dispersion Ub = 1 km/s for 
the background cluster. Here, the cross sections for binary 
star interactions are ~ 4.2 times larger than those for sin¬ 
gle stars at the high end of the eccentricity range e = 1. 
Similarly, the binary cross sections are ~ 3.2 times larger 
at the low end of the eccentricity range where e = 0.10. 
Averaged over the span of eccentricity considered here, 
the binary cross sections are larger by a factor of ~ 3.6. 
This factor decreases with increasing velocity dispersion, 
however, because the binary cross sections fall according 
to the relation (cr) oc , whereas the single star cross 
sections fall as (cr) oc With these scaling laws, the 

cross sections for binaries are only a factor of 2 larger 
(than those for single stars) when the velocity dispersion 
is increased to Ub ~ 20 km/s. 

These results can be interpreted as follows: At high 
asymptotic speeds, which occur for Ub ^ 20 km/s, the two 
members of a binary pass by the solar system quickly 
enough so that binary motion and planetary motion play 
only a minor role in the interaction (this speed is much 
larger than the mean orbital speed of either the binary 
or the outer planet). As a result, the two stars interact 
with the solar system in an almost independent manner, 
and the cross sections for binary interactions should be 
a factor of ~ 2 larger than those for single stars (for 
large Ub). On the other hand, lower impact speeds can 
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be comparable to the binary orbital speed and/or the 
planetary orbital speeds. In this regime, the motion of 
the binary stars relative to one another during the en¬ 
counter can increase their chances of interacting with the 
planets, thereby leading to larger cross sections. In ex¬ 
treme cases, resonant interactions can occur when the 
velocity scales of the prob lem are all comparable (see also 
iLaughlin fc AdamJ200ol l. and these long-lived events can 
greatly increase the chances of disruption of planetary 
orbits during the encounters. To be consistent with this 
picture, the ratio of the single-star cross section to the 
binary cross section must decrease less steeply with Vh, 
as found here. 

Before leaving this section, we briefly address the is¬ 
sue of how self-similarity can arise in the context of solar 
system scattering. In its full form, this problem has six 
velocities (four planetary orbits, one binary orbit, and the 
encounter velocity) and seven masses (four planets and 
three stars). One expects self-simi larity only when most 
of these scales do not contribute l|BarenbIattl l200,'j l . We 
can construct an argument to reduce the number of rele¬ 
vant scales as follows: To leading order — and only during 
the encounter itself — planetary interactions with the bi¬ 
nary are independent of interactions with other planets. 
As a result, we can (often) treat the encounters as single¬ 
planet systems scattering with binaries. The planet itself 
is usually small enough to be considered as a test mass, so 
that we are left with “only” three masses and three veloc¬ 
ities. The binary masses are always drawn from the same 
IMF, and the cross sections are determined through many 
samples of that IMF (Me ^ 80, 000), thereby leaving the 
ratio M*/(Mi* -|- M 2 *) as the most important mass vari¬ 
able. In the regime of interest, the cross sections have val¬ 
ues in the range (a) ~ 10'* —/ewx 10® AU^, which implies 
that the length scales that characterize the interactions 
Ic = ~ 100 — 500 AU. This size scale is larger than 

that of both the planetary orbits (a = 5 — 30 AU) and 
most binary orbits (where the peak of the period distribu¬ 
tions corresponds to as « 42 AU). If the orbital speeds of 
the planets and the binary are fast enough, then their or¬ 
bits can be replaced by rings of mass with the same sem i- 
major axis and eccentricity dMurrav fc Dermotd Il999l l. 
This averaging effectively eliminates the orbital veloci¬ 
ties from the problem and leaves the velocity dispersion 
Vb as the most important velocity variable. Indeed, we 
find that the cross sections depend most sensitively on 
the stellar host mass M* (equivalently, the mass ratio 
M*/(Mi* -I- M 2 *)) and the velocity dispersion Vb- This 
argument is not exact, however, and the additional scales 
(e.g., orbits speeds) do play some role. These complica¬ 
tions are responsible for the spread in the scaled cross 
sections shown in Figures [T^ - I16l 

We can also compare these scaling results to 
analytic results found in previous stud ies (see, e.g., 
iHeggie fc Rasid[i996l : ISnurzem et al.|[200^ , although the 
system parameters are not exactly the same. The lat¬ 
ter study finds a scaling relation {a) oc in the 

impulsive regime (where Uoo ~ Ub is much greater than 
the orbital speed of the planet) and (cr) oc for non- 
impulsive encounters. Our results (see Figures \T4\ and llTt 
are intermediate between these two scaling laws, since the 
encounters are rarely fully in the impulsive or the non- 


impulsive regime. In addition, this current study includes 
binaries, and the binary orbital speed is generally compa¬ 
rable to the planetary orbital speed. The binary motion 
can either add to or subtract from the relative velocity of 
the encounter (depending on the timing and geometry of 
the encounters), so that the scattering interactions have 
a wide range of relative velocities, even for a given Ub- As 
a result, our parameter space does not fall fully in any 
of the limiting regimes considered by previous analytic 
estimates. 


5 THE SOLAR BIRTH AGGREGATE 


Given that most stars are born within clusters, it is 
likely that the birth environment of our own Solar Sys¬ 
tem was a cluster of some type. The argument for a 
substantial birth cluster is bolstered by evidence for 
short-lived radionuclides in meteorites, which suggests 
that the e arly solar nebula was enr iched by a nearby 
supernova (Cameron fc TruranI 1 19771 : see the review of 
iDauphas fc Chaussidonll201lll . A number of previous pa¬ 
pers have considered how dynamical scattering encoun¬ 
ters in this putative birth cluster can provide constraints 
on the cluster properties (see the discussion of Section[T]). 
Unfortunately, however, no consensus has been reached. 
This section briefly revisits the issue in light of the up¬ 
dated cross sections determined above. 

The basic problem posed by the solar birth aggregate 
involves a number of ingredients: [I] Direct supernova en¬ 
richment of the early solar nebula requires a nearby mas¬ 
sive star, which is more likely to form in a larger stellar 
system. Further, significant nuclear enrichment requires 
close proximity (distances d = 0.1 — 0.3 pc), which im¬ 
plies that the supernova progenitor lives within the same 
cluster. Acting in the opposite direction, larger clusters 
can potentially disrupt planetary systems through the 
action of both [II] dynamical scattering (with the cross 
sections determined here) and through [III] intense radi¬ 
ation fields which can evaporate gaseous disks. In order 
for the solar system to reach its present-day state, how¬ 
ever, the orbits of the giant planets cannot be greatly 
perturbed and the early solar nebula could not be too 
severely evaporated. On the other hand, [IV] the classi¬ 
cal Kuiper belt has an apparent edge at ~ 50 AU, and 
[V] the dwarf planet Sedna has an unusual orbit; both 
of these solar system properties could be explained by 
requiring a close encounter with another member of the 
cluster. The challenge is to find a birth scenario for the 
solar system that successfully negotiates the compromises 
required to simultaneously explain all five of these con¬ 
straints. Supernova enrichment, the edge of the Kuiper 
belt, and the orbit of Sedna all argue in favor a large 
and long-lived cluster; disruption via both scattering and 
radiation argue in the opposite direction. 

Existing work has considered a variety of approaches 
to this issue. Several authors advocate solar birth clus¬ 
ters with stel lar membership size in the range N = 
10® --10* (e-K- - Adarns fc Lau'g hlinll200ll : |Portegies ZwartI 
I 2 OO 9 I : lAdamj l2010l : Pfalznerl 1201^ . These studies find 
that cluster systems in this decade of N lead to mod¬ 
erate dynamical disruption of their constituent planetary 
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systems. Additional work focuses on even larger, longer- 
lived clusters and find that they can instigate substan- 


tions (iMalmberg 

et al.ll2007l. 201ll; Smirzem et al.ll2009l: 

IParker & Ouanz 

2 OI 2 I; iHao et al.| |2013|l. On the other 


hand, competing work suggests that the solar birth clus¬ 
ter does n ot produce significant disruption of plane¬ 
tary orbits llWill iams fc Gaidosl l2007l : IPukes fc Krumhold 
I 2 OI 2 I : ICraig fc Krumholj|2013l : see also IWilliamsI 2010f) . 

The aforementioned papers thus reach different con¬ 
clusions about the importance of dynamical scatter¬ 
ing of planetary systems in clusters. These differences 
arise because of varying assumptions about cluster prop¬ 
erties and varying assumptions about how to enforce 
the five constraints on solar system properties outlined 
above. Althongh a full review of this topic is beyond the 
scope of this work, we provide a brief o verview below 
(for additional detail, see the reviews of Adamj I 2 OI 0 I : 
IPauphas fc Chaussidonll 201 ll : IPfalznerjl2ni3ll . 

For a given type of disruption, with cross section (a), 
the interaction rate is given by T = nt,{a)v (from equa¬ 
tion [ 1 ]). The total expected number A^dis of disruption 
events, per solar system, integrated over the lifetime r of 
the cluster is then given by 


[ rdt= [ 
Jo Jo 


Adis = / rdt= / nt{a)vdt. 


( 12 ) 


The number of disruptive interactions thus depends on 
the speed v at which a given solar system encounters 
passing binaries, their number density n*, and the total 
time T spent within the cluster. 

We first consider the speed v. Recall that the interac¬ 
tion cross section {a) varies with the velocity dispersion 
of the cluster according to the relation (a) oc {(7)ov^'^. 
If we identify the speed v with the velocity dispersion 
Ub of the cluster, then the product {ij)v oc . As 

a result, most of the velocity dependence of the cross 
section is compensated by that of the interaction rate, 
so that the number of disruption events depends only 
weakly on the velocity dispersion. As an example, con¬ 
sider the Orion Nebula Cluster (ONC), an intermediate¬ 
sized young stellar system with velocity dispersion Ub ~ 2 
km/s iHillenbrand fc Hartmannlll998h . Provided that it 
stays intact, the ONC is like ly to evolve into an open clus¬ 
ter resembling the Pleiades (iKrouna et aLllioOlh : over the 
coming ~ 100 Myr, the velocity dispersion of the cluster 
will slowly decrease to Ub ~ 1 km/s. Over this span of 
time, the quantity ' that dehnes the velocity depen¬ 
dence of the interaction rate varies by only about 32%. 

For setting the number of disruption events, one 
important quantity is the time r over which clusters 
remain intact as dynamical systems. In the simplest 
terms, although most stars are formed in clusters, these 
astronomical entities come in (at least) two distinctly 
different flavors. Only about 10 percent of the stel¬ 
lar population is born within clusters that are suffi- 
cently robust to become open clust ers llRobertsI 1 19571 : 
iBattinelli fc Canuzzo-Dolcettal Il99lh . which are rela¬ 
tively long-lived (r = 100 Myr - 1 Gyr). The remaining 90 
percent of the stel lar population is born within embed¬ 
ded clusters (e.g., I Allen et al.ll 200 ^ . which have much 
shorter lifetimes (r ~ 10 Myr). As shown below, solar 


systems that are born within long-lived clusters can have 
an appreciable chance of dynamical disruption; short¬ 
lived clusters lead to significant disruption with greatly 
reduced probability. 

Another important quantity is the density of the 
cluster. For clusters found in the solar neighborh ood, 
the cluster radius R oc (iLada fc Ladallioosh . so 

that the clusters displ ay nearly constant surface den¬ 
sity (I Adams et akllioOg l. With this relation, clusters with 
larger stellar membership sizes A have lower mean den¬ 
sities. However, the clusters in the sample are relatively 
small (with A < 2500), and this trend does not con- 
tinue up to the large st clusters with A = 10"^ — 10® 
ll Whitmore et akl 120071 1. or to the subpopulation of sys¬ 
tems that become globular clusters. The largest clusters 
can thus have larger densities. 

To assess the effects of scattering encounters, we need 
to specify the rate F at which solar systems encounter 

J iassing binaries (and single stars). As shown previously 
Adams et al.ll 200 ^ : IProszkow fc Adamj[2009ll . the rate 
of close encounters in a cluster can be written in the 
convenient form 

r = r,(l)h (13) 

where bo is a fiducial distance (taken here to be bo = 
1000 AU), and where the fiducial rate Fq and index 7 
depend on the cluster properties. The index 7 falls in the 
range 1 si 7 ^ 2 , where the extreme of the range corre¬ 
spond to perfect gravitational focusing (7 —>■ 1 ) and the 
full geometrical cross section (7 —>■ 2). In these systems, 
encounters beyond ~ 1000 AU are little affected by grav¬ 
itational focusing. Since the cross sections calculated in 
this paper include gravitational focusing, we can write 
the interaction rate in the form 

The benchmark interaction rate Fq has a typical value 
of about 0.1 interactions per target star per Myr. How¬ 
ever, given the wide range of possible cluster proper¬ 
ties, it can vary over a wide range, from an order of 
magnitude lower to an order of magnitude larger than 
this fi ducial value (see Figures 6 and 7, and Tables 8 - 
13 in IProszkow fc Adama 1200^ 1. Note that the bench¬ 
mark rate is, in general, larger than the simple estimate 
Fo ~ (n*)ub 7 r 5 o, where (n*) is the mean density of the 
cluster. The stellar density that defines the interaction 
rate is not the mean over the cluster, but rather the 
weighted mean over the integrated orbits of the ensemble 
of cluster members. The cluster members generally do 
not stay at a given cluster radius, and the cluster density 
is centrally concentrated, so that solar systems sample 
the higher stellar densities of the cluster core. This effect 
is amplihed by the starting conditions for clusters, which 
start with subvirial initial conditions; as a result, the or¬ 
bits are more radial than isotropic, resultin g in more ex¬ 
cursions through the dense cent ral core (see Adams et al.l 
I 2 OO 6 I : IProszkow fc Adam3l2009l for further discussion). 

Collecting the results outlined above, we can write 
the number of disruption events (from equation [T2]) in 
the form 
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The cross section for moderate solar system disruption 
can be taken as ((j)o ~ 2.5 x 10® AU^, which corre¬ 
sponds to events producing eccentricity increases Ae = 
0.1 and/or increases in the spread of inclination angles 
Ai = 10° (e.g., see Figures [T] and O. To obtain this 
value, we use a linear combination of the binary and 
single-star cross section (see Figure ITlTl . and an assumed 
binary faction of 2/3. Although these changes to the or¬ 
bital elements are not devastating, they are large enough 
to distinguish a disrupted solar system from our own. 
Note that this value can be written ((t)o ~ (500 AU)^, 
which is somewhat larger than the previous estimate of 
~ (400 AU)^ Ifrom lAdams &: LaughliniboOll lFI The lead¬ 
ing factor in equation (HU is thus of order 1/10. Since the 
benchmark interaction rate Fq ~ 0.1 Myr”*^, the cluster 
lifetime must be relatively long, r ~ 100 Myr, in order for 
disruption to take place with high probability. In other 
words, most solar systems residing in long-lived clusters 
can experience moderate disruption. 

The cross sections for planet ejection are smaller 
than the values used above by a factor of ~ 3. As a result, 
only a fraction (~ 1/3) of the solar systems in long-lived 
clusters are expected to lose planets with Neptune-like or¬ 
bits (with even smaller fractions for closer planets). Keep 
in mind, however, that the benchmark interaction rates 
Fo can vary by a factor of ~ 10 in both directions. 

The above considerations resolve some of the dif¬ 
ferences found in the literature concerning the disrup¬ 
tion rates for planetary systems in clusters. In order for 
disruption to occur with high probability, clusters must 
live for relatively long times r ^ 100 Myr. Indeed, the 
studies that find low disruption rates consider the clus- 
ters to have relatively short lifetimes r ~ 10 Myr (e.g., 
IWilliams fc Gaid^l2007l : iDukes fc Krumh^l2012h . 

How long are clusters expected to stay together? 
As outlined above, the cluster population has at least 
two branches. Some clusters disperse over relatively short 
time scales of only ~ 10 Myr. The robust clusters that 
survive to become open clusters have empirically deter¬ 
mined lifetimes that can fit with a function of the 
form 

= 2.3Myr j , (16) 

where Me is the cluster mass dLamers et al. I l2005ll . With 
this relation, clusters with initial masses larger than 
~ 550Mq live longer than 100 Myr and can potentially 
disrupt their constituent solar systems. More specicifally, 
we can write the dynamical constraint in the form 

fVdis « ■^(ro)2.3(TiV)®/®((ub/lkm/s)-=/®)^l, (17) 
ttOq 

^ The difference arises because the present study increases the 
target area in equation (O from B = i to B = 100, thereby 
including more distant e vents. Note that the original work 
ll Adams &: Laughllnll200lh correctly introduced the cross sec¬ 
tions as lower limits. The present cross sections are also lower 
limits, although they are much closer to their greatest lower 
bounds. 


where T (~ 1/2 Mq/ stax) is the mass-to-number ratio 
(the conversion factor between cluster mass Me and clus¬ 
ter membership size N), and where we include the time 
average of the velocity dispersion of the cluster (raised 
to the proper power). After some rearrangement and the 
specification of typical numbers, this constraint can be 
written in the form 


N < 5000 


(^)o ^ [ (ro) \ 

2.5 X 105AUV V0.05Myr-V 


(18) 


Vb 


1 km/s 


- 2 / 5 \ 1 - 5/3 


^ 10 ^ 


Note that the disruption cross section is determined more 
precisely than either the expected age of the cluster (from 
equation m) or the benchmark interaction rate Fq. This 
latter quantity can be determined to high accuracy for 
a given set of cluster properties and initial conditions, 
but its value varies appre ciably from cluster to cluster 
JProszkow fc Adams! 12009 ! '). Equation (1181) uses a value 
near the low end of the range in order to provide an 
upper limit on N. In light of these uncertainties, a rea¬ 
sonable order-of-magntidue estimate for the dynamical 
constraint is Al ^ 10"^, as given by the final inequality. 
This result is roughly con s istent with previous e s timates 
jAdarns fc Laug hlinlpOOll : IPorteeies Zwart|[2009l : lAdamj 
120101 : Pfalznei 20131) . I^netheless. the full probability 
distribution for the survival (or disruption) of planetary 
systems as a function of cluster size N should be con¬ 
structed. 

The constraint given by equation nil) assumes that 
the solar birth cluster is relatively long-lived. If the so¬ 
lar system formed within a cluster that disperses in only 
~ 10 Myr, the corresponding dynamical constraint would 
be considerably weaker. The motivation for considering 
a long-lived cluster comes from constraints jointly im¬ 
plied by the five solar system properties outlined at the 
beginning of this section. Direct supernova enrichment 
[I] favors a long-lived cluster, so that the progenitor star 
has enough time to live, evolve, and explode. An even 
stonger argument comes from the need for a scattering 
event to produce the edge of the classical Kuiper belt [IV] 
and to produce the orbit of Sedna [V]. If these solar sys¬ 
tem properties arise from dynamical interactions in the 
birth cluster, then a long-lived stellar system is strongly 
indicated. It remains possible for these features of the 
solar system to be explained in other ways. Nonetheless, 
any self-consistent set of constraints on the solar birth 
environment must explain all three of these properties, 
and must simultaneously account for the corresponding 
constraints due to dynamical s cattering encounters [II] 
and radiation fie lds [III] (e.g., see lFatuzzo fc Adamsl2008l : 
lThomps^l2013l~l . 

For completeness, we also consider possible con¬ 
straints on the solar birth cluster for the scenario where 
the solar system spends much of its early life in the ultra¬ 
compact multi-resonant configuration (see Section|3] Fig¬ 
ure ID). The cross section for removing the solar system 
from its resonant state is then given by equation m, 
which is more than nine times larger than that used 
above. If, in addition, the removal of the solar system 
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from resonance always led to significant disruption over 
longer times, then the maximum size of the solar birth 
cluster would be ~ 40 times smaller than that of equa¬ 
tion In practice, however, the solar system, after 

being removed from resonance, will not always be sig¬ 
nificantly disrupted (for example by ejecting a planet) 
before it ev olves and spreads out (as advocated by the 
Nice model: lOomes et'Sl200,‘il : iTsiganis et al.ll200,’ih . To 
assess the risk of disruption in this case, one must also 
know the probability of the non-resonant (but still com¬ 
pact) solar system experiencing disruption on sufficiently 
short time scales. This calculation involves a large ensem¬ 
ble of long-term (~ 100 Myr) solar system integrations 
and is beyond the scope of this present work. Nonethe¬ 
less, direct application of equations (0) and suggests 
that the constraint could be tighter than that derived for 
the solar system in its usual configuration. 

Finally, we note that another class of observational 
constraints on the solar birth environment might become 
available. Given that the birth cluster is expected to have 
N « 10® —10"^ stars with similar chemical composition, it 
is possible in principle to find other members of our solar 
birth aggregate. Although billions of years have passed 
and the cluster has long since dispersed, perhaps ~ 20 
of th ese solar siblings could reside within 100 pc of the 
Sun (IPortegies Zward 120091 '). By focusing on the chemi¬ 
cal species that show the most variation from cluster to 
cluster, it is possible to o bservationally disting uish these 
siblings from other stars llRamfrez et al.l[2014l i. The dis¬ 
covery of even a few such stars would provide strong con¬ 
straints on the properties of the solar birth cluster and its 
location within the Galaxy. On the other hand, the Solar 
System could have had a more complicated dynamical 
history including large radial migration in the Galaxy 
dKaib et al.l [20111 1. which could reduce the chances of 
finding solar siblings. 


6 CONCLUSION 

6.1 Summary of Results 

Using results from more than 2 million individual numer¬ 
ical scattering experiments, this paper has found cross 
sections for the disruption of planetary orbits in solar 
systems interacting with passing stars and binaries. Our 
specific results can be summarized as follows: 

[1] More compact solar systems have smaller inter¬ 
action cross sections (Figures [T] and [2|. To leading order, 
the cross section for a given disruption event (e.g., planet 
ejection or eccentricity increase) scales linearly with the 
semimajor axis of the initial orbit, i.e., (a) oc a (see Fig¬ 
ure [TH). 

[2] For most solar systems, the cross section for a 
given planetary orbit to be disrupted during a scatter¬ 
ing encounter is almost independent of the other planets. 
This feature of the interactions allows for the scaling anal¬ 
ysis presented in Section|T] Of course, after the encounter, 
solar systems that suffer moderate disruption can subse¬ 
quently experience orbital instability, and this latter ef¬ 
fect does depend (quite sensitively) on the other planets 
in the system. In addition, for highly self-interacting solar 


systems, those with sufficiently massive planets and/or 
close orbits, interactions among the planets themselves 
can lead to effectively higher cross sections (e.g., see Fig¬ 
ures [5] and 0). 

[3] The dependence of the cross sections (cr) on the 
post-encounter eccentricity e has a nearly exponential 
form (see Figures [1]-[5}. As a result, the cross sections 
can be written (cr) oc exp[—fee], where fe ~ 4/3 provides a 
good fit across the range of parameter space considered 
in this work (Figures fTHlfTTl andlT^. 

[4] The cross sections depend sensitively on the veloc¬ 
ity dispersion Vh of the background environment, where 
the dependence displays a nearly power-law form. More¬ 
over, the shape of the cross section curves, as a function 
of eccentricity, are nearly the same across the parameter 
space considered here (Figures 0 and [TT]|. The cross sec¬ 
tions can thus be written as (cr) oc exp[—fee], where 
7 = 7/5 and fe = 4/3 provide a good fit over range of 
interest. 

[5] The cross sections depend on the mass M* of the 

host star, where the dependence has the approximate 
form (cr) oc The mass dependence is somewhat 

more complicated, however, as the cross sections are not 
fully self-similar (see Figures [4] and [15}. For more disrup¬ 
tive encounters (where e 1 and planets are ejected), 
the scaling with mass is somewhat steeper and the form 
(cr) oc pr ovides a better fit (c onsistent with previ¬ 

ous results from [Adams et ^l2006ll . 

[ 6 ] Most of this work considers planetary orbits with 
vanishing initial eccentricity e. Nonetheless, for solar sys¬ 
tems starting with e 7 ^ 0 , the interaction cross sections for 
eccentricity increase are nearly the same (Figure 0, pro¬ 
vided that one considers post-encounter eccentricities suf¬ 
ficiently larger than the starting values (roughly, by the 
increment fee ~ 0.1). This finding stands in contrast to 
the related problem of single stars interacting with bina¬ 
ries, where the cross sections for binaries with e = 0 a nd 
e 7 ^ 0 are significantly different (iHeggie fc Rasiol[l996l i. 

[7] The above results can be combined to write the 
cross section for eccentricity increase for solar systems 
interacting with binaries in the general form 


(a) = 4050 (AU)^ 


( “ 

(M.\ 

VAU/ 

\MeJ 


-1/3 


(19) 


X 


Ikm/s 


-7/5 

exp 




This result holds over the ranges of parameters given by 
5 AU < a s; 50 AU, O.25M0 ^ M* sj 2Mq, 1 km/s 
^ Ub ^ 16 km/s, and 0.1 ^ e ^ 1. Equation (1191) is in 
good agreement with the numerically obtained results: 
For fixed stellar mass, the RMS relative error for the 
range of starting semimajor axis, velocity dispersion, and 
post-encounter eccentricity is less than about 12 percent 
(see Figures [ 13 ] and [14}. Including variations in the stellar 
mass, the RMS error is less than about 20 percent (Figure 
IS}. Over the same regime of parameter space, the cross 
section itself varies by more than a factor of ~ 1000. 
Equation IS} provides the total ejection cross sections 
(including capture events) in the limit e —>■ 1; the cross 
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sections for ejection and capture are listed separately in 
Table 1. 

[8] The cross sections for increasing the spread of 

inclination angles Ai are comparable to those for in¬ 
creasing eccentricity (Figure [7|). The cross sections for 
Ai also show a nearly self-similar form, and scale with 
velocity dispersion of the background cluster according 
to (cr) oc fFigure fTUll . This scaling exponent is the 

same as that found for eccentricity increases. The cross 
sections can be fit with an exponential dependence on the 
variable x = sinAz. Although inexact, one can identify 
increases in inclination with increases in eccentricity such 
that Ax ~ Ae. In general, increases in the spread of incli¬ 
nation angles and orbital eccentricity are well-correlated 
(Figure[8}, although the Ai values for a given Ae display 
a wide range. We have also determined the cross sections 
for increasing the inclination angles beyond 39.2°, the 
benchmark value required for the Kozai effect to operate 
(equation [TO]'). 

[9] In addition to the ejection of planets during the 
scattering encounters, orbital eccentricites can be in¬ 
creased so that planetary orbits will cross each other. 
Most solar systems in such states will eject — or perhaps 
accrete — planets on relatively short time scales. For sys¬ 
tems with the architecture of the current solar system, 
the cross section for this channel of secondary ejection is 
comparable to that of direction ejection, so that the total 
cross section for ejection is effectively doubled (Figure^. 
For the ultra-compact configuration of the solar system 
(in or near multiple mean motion resonances), the cross 
section for ejection due to orbit crossing is comparable to 
that of the standard solar system, but the cross section 
for direct ejection is smaller. 

[10] The cross sections for changing the semimajor 
axes of the planetary orbits are smaller than those for 
increasing eccentricity and/or inclination angle (Figure 
[TOJ. Equivalently, the semimajor axes change much less 
than the other orbital elements during scattering encoun¬ 
ters. In rough terms, 10% changes in the semimajor axis 
— for planets that remain bound — have approximately 
the same cross sections as planetary ejection. 

[11] The cross sections for solar systems interacting 
with single stars are smaller than those for binary encoun¬ 
ters [Figure [TTIl . The single-star cross sections are nearly 
self-similar (Figure fTTll . and scale with the semimajor axis 
of the planet and cluster velocity dispersion according to 
(a) oc a The scaling exponent for velocity is some¬ 
what smaller than that for binaries and the dependence 
of the cross sections on the post-encounter eccentricity 
is steeper. On average, the single-star cross sections are 
smaller than the binary cross sections by a factor of ~ 3.6 
for small velocity dispersions (ub = 1 km/s). This factor 
falls to only ~ 2 for larger values Vh ~ 20 km/s; for 
higher speeds we expect the binary components to act as 
two separate stars during the encounters (except for close 
binaries). In general, the effective cross section is a linear 
combination of the single and binary star cross sections, 

O ~ fbOhiu3.Ty T (1 /b) (^)single 5 (^0) 

where fb is the binary fraction. 

[12] We have briefly revisited the dynamical con¬ 
straint that can be placed on the birth aggregate of our 


solar system due to scattering encounters (Section[5]). The 
strength of this constraint depends crucially on whether 
one assumes that the solar system forms in a robust, long- 
lived cluster (with r > 100 Myr, like those that become 
open clusters) or in a short-lived cluster that dissipates 
within r ~ 10 Myr. For long-lived clusters, the require¬ 
ment that the solar system is not disrupted implies an 
order of magnitude upper limit on the solar birth ag¬ 
gregate of A^ 10^ (see equation [18]). In practice, one 
should construct the probability distribution for solar 
system survival/disruption as a function of N (using the 
cross sections determined herein), and combine it with 
the other con straints on the birth cluster (se e Fi gure 7 in 
Adams!l2010l : see also jPortegies Zwartll20()^ and lPfalznerl 
20131'! . 

[13] The cross section for removing a solar system 
from mean motion resonance is much higher than that 
required to disrupt the planetary orbits. For the ultra¬ 
compact multi-resonant conhguration advocated by some 
versions of the Nice model, this cross section (see equa¬ 
tion El) is ~ 9 times larger than the disruption cross 
section for the usual solar system architecture. If removal 
from resonance leads to longer-term instability, then con¬ 
straints on the solar birth aggregate would be tighter for 
systems in the multi-resonant configuration. 


6.2 Discussion 

The cross sections reported in this paper are subject to 
three different types of uncertainties, and the distinctions 
among these quantities should be kept in mind. [1] First, 
the Monte Carlo procedure used to determine specific 
cross sections (as outlined in Section O results in uncer¬ 
tainties due to incomplete sampling. These uncertainties 
decrease with increasing size of the ensemble of simula¬ 
tions and are proportional to . Over most of the 

parameter space, we run sufficient numbers Me of scat¬ 
tering experiments so that the sampling errors are less 
than ~ 5% and usually even smaller. These sampling 
errors are present in all of the cross sections presented 
in Section |31 although they are usually not included on 
the plots (however, see Figure [5]|. [2] Next, in Section (4] 
we explore scaling laws to collapse the cross sections for 
varying velocity dispersion vt, host mass M*, and planet 
semimajor axis a into nearly self-similar forms. The range 
of the resulting scaled functions is thus characterized by 
the error bars shown in Figures 1131114111511161 and 1171 
These error bars represent a measure of the degree to 
which the cross sections depart from self-similarity. The 
size of these error bars falls in the range 10 - 15%, ex¬ 
cept for the scaling with the mass of the host star (where 
the error bars correspond to 20% departures). [3] Finally, 
the mean of the scaled cross sections are described by ht- 
ted functions with simple forms. The differences between 
these functions and the mean scaled cross sections are of 
order 5 - 10%, smaller than the standard deviations of 
the different sets of cross sections used to construct the 
mean forms. 

In addition to the uncertainties outlined above, the 
cross sections calculated herein depend on the features of 
the stellar population that provides the perturbations. As 
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described in Section the cross sections sample the dis¬ 
tributions of stellar masses, binary periods, binary mass 
ratios, binary orbital eccentricities, etc. Different choices 
for these distributions will lead to corresponding varia¬ 
tions in the cross sections. Although we use observations 
to specify the distributions, they are nonetheless subject 
to both measurement error and possible variations from 
region to region. 

The numerical simulations carried out for this paper 
determine the immediate changes in the orbital elements 
of the solar systems due to passing stars. However, ad¬ 
ditional changes in the orbital elements can occur over 
longer time scales. As one example, after an encounter, 
a planetary system often has larger eccentricities, which 
can lead to orbital instability over longer spans of time. 
But the timescales for such instabilities can have a wide 
range. For systems where the eccentricities are increased 
so much that planetary orbits cross, one expects insta¬ 
bility and (usually) planet ejection on a relatively short 
time. The cross sections for orbit crossing are thus of 
great interest and are given in Figure[9] For systems with 
smaller eccenticity increases, however, orbital instability 
can take much longer. For compact multi-resonant solar 
systems, modest changes in the orbital elements and/or 
the removal of the system from its reso nant state can lead 
to in stabilities over millions of years dSatygin 8z Bro\^ 
I 2 OIOI') . For systems with more widely separated orbits, in¬ 
stabilities can_take^yenJongCTthM the current age of the 
unive rse llBatvein fc LaughliDll20()§ : iLaskar fc Gastineaul 
I 2 OO 9 I '). To study this issue, the post-encounter solar sys¬ 
tems must be integrated over long time scales (up to 
billions of years) to fully determine the effects of the 
encounters. This task is left for future work. On an¬ 
other front, the orbits could also damped after the scat¬ 
tering encounters, thereby moving the orbits back to¬ 
wards smaller eccentricit ies dLevison fc MorbidellillioO^ : 
IPicQgna fc Marzarill201 ^ 1. This effect should also be con¬ 
sidered in follow-up studies, especially on time scales of 
1 - 100 Myr when solar systems are expected to retain a 
significant population of planetesimals. 

The scattering encounters considered herein can be 
effective in sculpting giant pl anet orbits and the Kuipe r 
Belt of our Solar System fe.g.. lKenvon fc Bromlevir2nn4l 'l . 
On the other hand, the Oort cloud is too large to be 
produced within a young embedded cluster (e.g., see 
iBrasser et al.l 12012 for further discussion). More s pecif¬ 
ically , the Oort clo ud extends out to ~ 50, 000 AU (lOortI 
Il95nl:l,lewidl2nnil'l . more than 1000 times the size of the 
solar systems considered in this paper. With this enor¬ 
mous size, the Oort cloud would be decimated by pass¬ 
ing stars within the cluster. As a result, the cloud must 
be produced later, after the solar system leaves its birth 
cluster, or perhaps during its exit. Any viable scenario 
for the solar birth environment must simultaneously ac¬ 
count for the Oort cloud, the giant planet orbits, Kuiper 
Belt properties, radioactive enrichment, Sedna’s orbit, 
and survival of the solar nebula gas reservoir; these cou¬ 
pled constraints thus pose an interesting and challenging 
opimization problem for further study. 
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